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Abstract—To design rationally a chemical reactor with its associated control equipment, it is 
necessary to study the transient behaviour of the system. In this series of papers a well-agitated 
continuous reactor either of the liquid or fluidized variety is considered under various types 
of control and criteria are derived enabling one to determine under what conditions different 
modes of control are successful. This first part, after an introduction to the whole subject, 
introduces the method of linear approximation and goes on to discuss the possibility of con- 
trolling a naturally unstable steady state by various systems of perfect and imperfect control. 


Résumé—L’étude rationnelle d'un réacteur chimique avec son appareillage de contrdle, 
nécessite la connaissance du comportement transitoire du systéme. Dans cette série d’articles, 
les auteurs considérent avec différents types de contréle, un réacteur continu bien agité, soit a 
liquide, fluide ou variétés fluidisées. On obtient ainsi des critéres permettant de déterminer dans 
quelles conditions ces divers modes de contréle conviennent le mieux. Aprés une introduction 
relative 4 l'ensemble du sujet, cette premiére partie présente une méthode d’approximation linéaire, 
et se poursuit par la discussion de la possibilité de contréler un état de régime instable naturelle- 
ment, par divers systémes de contréle. 


Zusammenfassung— Zur rationellen Berechnung eines chemischen Reaktors und der zuge- 
hérigen Regeleinrichtung ist es notwendig, das nicht-stationiire Verhalten des Systems zu 
untersuchen. In dieser Reihe von Veréffentlichungen soll ein gut geriihrter kontinuierlicher 
Reaktor mit fliissiger Phase oder als Fliessbett unter verschiedenen Regelungsarten betrachtet 
werden. Dabei will man die Kriterien kennenlernen, nach denen es sich entscheiden lisst, unter 
welchen Bedingungen die verschiedenen Regelungsarten erfolgreich arbeiten. Dieser erste Teil 
gibt, nach einer Einfihrung in den gesamten Gegenstand, einen Cberblick tiber die Methode der 
linearen Approximation und diskutiert dann die Méglichkeit, einen von Natur aus unstabilen 
stationiren Zustand durch verschiedene Systeme volstindiger und unvollstindiger Regelung 
zu regeln. 


1. INTRODUCTION all of these are stable. The first treatment of the 


Tue behaviour of the continuous flow stirred tank 
reactor has received considerable attention since 
it was first discussed by MacMULLIN and WEBER 
[9] and more particularly since DenBicn’s paper 
in 1944 [5]. Van Heerpen [12] showed how 
certain steady states are possible at which the 
reactor behaves autothermically and that not 


transient behaviour, taking into account the 
essential non-linearity of the problem, was given 
in 1955, [1] and [2]. The methods of non-linear 
mechanics, which have been developed over the 
past 50 years, were used to study the approach 
to the steady states and their stability. It is the 
purpose here to use these methods to study the 
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problem of controlling such a reactor at a steady 
state. They prove admirably suited to the 
discussion of this problem and some care is taken 
to present them in a manner appealing to the 
chemical engineer. 

The two principal notions involved are those 
of the locally linearized equations (par. 3) and the 
phase plane (par. 8). The former is a set of linear 
equations which approximate to the non-linear 
ones in the immediate vicinity of a steady state 
and allow the question of its stability to be 
determined. The latter is the representation 
of the behaviour of the reactor as a curve in a 
plane with co-ordinates temperature and com- 
position and this shows if and how the steady 
state is approached from given initial conditions. 
The former notion is used in par. 5 to decide what 
form of control will be successful and in par. 9 to 
show the nature of the steady states during the 
The latter, 
arguments 


evolution of a successful control. 


together with certain topological 
based on it, is used in par. 9 to show the complex 
nature of the control problem and to predict 
the phase plane configuration. In par. 10 another 
tool of non-linear mechanics is introduced, the 
so-called bifureation theory of Porxncare [11] 
to discuss the interesting phenomenon of limited 
control which may arise even when the controlled 
These 


exemplified and confirmed by an extensive series 


steady state is stable. methods are 


of calculations described in par. 11. Further 
topics necessary to tie in the loose ends of this 
discussion are presented in the remaining para- 
graphs : p. 5 imperfect control, such as arises 
by the time lags in an actual servo-mechanism ; 
par. 13, off-on control, the limit at very high 
amplifications when an instrument meets the 
limits of its movement almost instantaneously; 
par. 14, the extension of the discussion to many 
reactions taking place simultaneously in a series 


of reactors. 


2. Tue EQvuATIONS 


The main concern will be with a single con-’ 


tinuous agitated tank reactor. This is a vessel 


of volume V into which a stream of reactants 
llows at a constant flow rate g. The vessel is 


stirred so that the composition and temperature 


of its contents are constant throughout and these 
are the composition and temperature of the 
product stream. For the moment attention will 
be confined to a single reaction taking place 
between n chemical species and involving no 
Such a 


change of volume. reaction 


written 


may be 


0 (2.1) 


where one of the reactants to be designated by X 
without a subscript has been singled out. This 
choice is arbitrary, but if it is chosen to be that 
reactant which would be first exhausted if the 
reaction were not limited by equilibrium then 
the concentration # of this substance cannot be 
Let 
av, and a, be the concentrations (in moles per unit 
the 
species X,, T be the temperature of the reactor 


greater than its influent concentration 2). 


volume) and influent concentrations of 


and ¢ the time. Then a unique reaction rate R 
can be defined, which will in general be a function 


of all concentrations and the temperature, by 


1 /dzx, d 
4, \a), , (i) 


r 


d 
where | ») denotes the rate of change of concen- 
( r 
tration due to reaction only. This reaction rate is 
unique for by the stoicheiometry of the reaction 


(2.2) 


(g=1,2,...% 1) 


1 dx, da 
dt 


so that R can be uniquely defined as a reaction 
rate in moles per unit volume per unit time. 
A mass balance for each substance gives the 
equations 
da 


V a)— VR (2.3) 


q (%» — 2) VAR (2.3a) 


but, by substituting from (2.2) in (2.8a), (2.3) is 
obtained so that this is the only equation needed, 
Moreover one can substitute in R (2;,2, T) the 
(2.2 have 


expressions for a; given by and 
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R= R(x, T; %, %) = R(x, t) in which ay and 


rg appear only as parameters. 


Let (— AH) = heat evolved mole of X 


formed 


per 


specific heat of reacting mixture 
density of reacting mixture 
temperature of reacting mixture 
inlet temperature of reactants 
heat 
from the reactor by the cooling 
U* is a function of T 


and the temperature and flow 


per unit volume removed 


svstem. 


rate of the coolant. 


Then a heat balance for the reactor gives 
_ Ox 
leo 
it 


( 


gep (T,,—T) — VU*—(—AH)VR (2.4) 


To avoid undue complication it is assumed that 


the thermal properties and inlet composition of 


the reactants are sensibly constant. This is no 


essential restriction and it could be removed 


at any time, but it does save carrying along a 
number of rather slowly varying functions of 
temperature. 

unit of time is 


A natural og the so-called 


holding time of the reactor; similarly x, is a 
useful unit of concentration and if Y is chosen as 
suggested above (a ZX») never exceeds unity; the 
natural unit for temperature is the difference in 


AH Ly Cp, 


from the complete conversion of YX, 


result 
Thus the 


temperature, which would 


following substitutions 


T qt reduced time 
€ r/Zo reduced concentration 
cpT ( AH) x, 
P VR gry 
Uv VU* qa, (— All) 


reduced temperature 
reduced reaction rate 


reduced heat removal 
reduce the mass and energy balance to 


lé ; . 
— =i - é P (é, ») 
dr 


(2.3) 


dy 


7. = (no — 9) — P(E) —- Ulm) (2-6) 
aT 


the 
concentration and temperature are such as to 


The steady state values £, and », of 


12: 


These 
may be found by first of all solving € = 0 for 


make £ and 7 vanish simultaneously. 


£ as a function of 7, say » = £,(7), and then 
substituting in the equation » = 0. Re-arranging, 
there results 


- U (n) P (€,(n), ») (2.7) 


7) No 


This equation may be given a physical interpre- 
tation as is done by Van Heerpen [12]. For 
writing it in the original variables 
qgep (T — Ty) Vu* (— AH) VR 

the left hand side is the total rate of heat removal 
both in warming up the reactants and by the 
cooling system. The right hand side gives the 
rate of heat generation and for a steady state 
the two must be equal. In general, the left hand 
side of equation (2.7) does not depart much from 
a linear function of » such as is shown by the 
line L in Fig. 1. 


sigmoid in shape (ef. the curve R in Fig. 1) 


The right hand side is generally 


though it has been shown [2] that more complex 
shapes are possible. Thus in cases commonly 
encountered there are three steady states corres- 
ponding to the intersections A, B, C of L and R 


in Fig. 1. 





. 1. Steady states with and without control. 


The steady state corresponding to the point B 
in Fig. 1 is clearly unstable, for, for slightly higher 
temperatures the heat generated is greater than 
the heat 
continue to rise until the steady state C 


absorbed and temperature would 
were 
reached. Similarly a slight decrease of tempera- 
ture would send the system to the steady state 4. 


It is evidently necessary for stability that the 
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gradient of L should be greater than that of R 
at the point of intersection. 
From (2.5) 


(1 + P,) dé, +P =0 


ly ’ 


(the notation P, = is being used) so that 


d DE 2) - - p, abs _ P=-—P (i 
dy dy : 


It is therefore necessary for stability that 


(1 Bi Pj t U,) ‘ P. 


P,). 


(2.8) 


It will be shown in the next paragraph that this 
condition is not sufficient, it is further necessary 
for stability that 2+ P,+ U,>—P, but 
this latter condition does not seem to have an 
obvious physical meaning. 


8. LocaL LINEARIZATION 


It was shown by Liapounorr [8] that questions 
of stability of a steady state could be determined 
by considering only the linear equations describ- 
ing small departures from the steady state. 
This theorem has been applied to the stirred 
tank reactor [1, 2] and the argument presented 
here simply applies this to the more general 
system under discussion. 

The equations (2.5) and (2.6) are of the form 


g P(é,7), 7 = Q (£, ») 


and the steady state £7, is 
P (é,, n,) = Q(é,. 7.) = 0. In the 


hood of this steady state one may write 


(3.1) 


that 
neighbour- 


such 


d 
ai (€ — é,) 
(P,), (€ — €,) + (P,), (4 — 0.) + 


7 () - 5) ai 


d 
d 





(Q.), (€ i €,) T ()s (Qy = 15) , 0, 
P 


where (P,), is - 


and O, 
d€ 


evaluated at £,, », ete. 


denotes terms of the order of squares and products 
of (€ — €,) and (» — ,). The linear system of 
equations obtained by neglecting the O, terms 
is stable only if 


— (P,), — (Q,), > 9 (3.3) 


(P.), (Q,), (3.4) 


To apply these criteria, three dimensionless 
numbers which characterize the steady state are 


(P.), (Qe), >o0 


introduced. They are : 
V aR 
q Mw 
increase in reaction rate with increasing z 


= (P,), 


increase in flow rate of 2 with increasing & 
. V oU* 
B=(U,), = 
qcp v7 
increase of heat rejection rate per unit temp. rise 
increase of heat content flow per unit temp. rise 
V(—AH) oR 
qcp MT 


increase of heat envolved by reaction per unit temp. rise 


(— P), 


increase of heat content flow per unit temp. rise 


Then applying the conditions (3.3) and (3.4) 
to equations (2.5) and (2.6) stability requires 


l+a+1+f>y (3.5) 
(l+a)(l+8)>¥y (3.6) 


These conditions may be given a geometrical 
representation if a, 8, y are thought of as co- 
ordinates of a point S in x, 8, y space. The surface 


E detined by 


is a plane intersecting the co-ordinate planes 
x = Oand £ = 0 in the straight lines y = 2 + 8, 
y = 2+ arespectively. The surface F defined by 


d,=(1 


/ x(1+p)—-y=0 (8.8) 


is an hyperboloid intersecting planes of constant 
x or 8 in straight lines. The quantities d, and d, 
are the distances from the point S to the surfaces 
E and F respectively and the condition for 
stability is that these should both be positive. 
It will be convenient to divide the «, 8, y space 
up into one stable and three unstable regions as 
follows 


Region 





Sign of d, 
Sign of d, 
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More information about the behaviour in the neighbour- 
hood of a steady state can be extracted from this local 
linearization. The solution of equations (3.2) is 


'—¢& 
7—% 


Cc; exp wy t C3 exp we t (3.9) 


Cy exp at Cy exp wet 


where w, and wg are roots of the quadratic 


° 
- 


w d, w dy 0 (3.10) 


and the constants Cy and Cy are not independent of C, 
and Cg. The stability condition d, > 0, dy > 0 simply 
ensures that the real parts of @, and wg shall be negative. 
If w, and ws are complex conjugates the right hand sides 
of equations (3.9) are damped oscillations and might be 
written in the form 


£ f Ay sin (pl ¢,) exp (AD) 


” le 1. sin (jal ¢g) exp (Al) 


When the approach to the steady state is of this oscillatory 
character the critical point £ 
it is stable or unstable according as A - 


7, is called a focus (Fig. 2a); 
Oo or A> oO, 


The stable case corresponds exactly to the motion of a 
simple harmonic pendulum with less than critical damping; 
it performs an infinite number of oscillations but their 
amplitude decreases exponentially with time (Fig. 2b). 
If A= 0 (ie., d, 
origin would persist as a simple harmonic oscillation, 


0) a small disturbance around the 


neither tending to grow nor decreasing in amplitude. The 
critical point is then called a vortex, and it stands at a 
If both roots w 
are real and of the same sign, the motion has no oscillatory 
character and corresponds to an overdamped pendulum. 


transition from stability to instability. 


It is stable or unstable according as the roots are negative 
or positive. The critical point is known as a node in this 
case (Fig. 2c, d). 








ts 7] 


Fic. 2. Types of critical point. 


If both roots are real but of opposite sign the critical 
point is called a saddle-point or col. Making use of the 
dependence of Cg and Cy on C, and Cg alluded to above, 
the solution of the linearized equations may be written 


&—- 
u] "s 
(l+«a 


Cy y exp ( 
C, (1 - @ 


w, b) Cg y exp (+ we t) 


w,) exp ( ~ wy t) 4 Ce 
wg) Exp (+ we t) 


where now w, and wg are both positive. Thus if originally 
(é £,)/(n — 5) y/(l+« w,) one may set Cy = 0 
and both £ and » tend exponentially towards £ and », 
and the ratio (€ — €,)/(m — ,) does not change (curves 
A and B of Fig. 2e). On the other hand if (£ — £,)/(m — n,) 
wg) originally £ and » tend exponentially 
away from £,, n, (curves C and D of Fig. 2e). 


y/l+a 
From any 
other original point Cg cannot be zero so that though 
&— &, j, may decrease at first the positive ex- 
ponential ultimately dominates and the condition departs 
from the steady state. In the ¢, » plane the trajectories 
appear as in Fig. 2e and the curves 4B and CD are known 


or 7 


as the separatrices. 
The condition that both the roots of equation (3.10) 
should be real is that d? 


~ 4 dy i.e. 


By? —2(a Bby+7>O 


(3.11) 


° ‘ » 
a* + p2 + > 


»2(a8 + By + ya) 


In the «, 8, y space this is the exterior of a parabolic 
surface touching the planes x« = 0 and y 
8 and y = « and the surface F along its inter- 


section with the plane EZ, 


0 along the 
lines y 
A section of these surfaces by a 
plane of constant 8 is shown in Fig. 3. 








Fic. 3. Nature of the critical points as functions of the 


parameters «, 8, y. 


4. Mopes or PERFECT CONTROL 
The differential equations (2.5) and (2.6) 


contain two variables 


€, the concentration 
n, the temperature 
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and a number of parameters 


£io9, the inlet concentrations 

MH, the inlet temperature 

y.» the temperature of the coolant 
q., the flow rate of the coolant, 


the last two appearing in the function U (»). 
At a given steady state all these quantities have 
fixed values and it is proposed to control this 
steady state by measuring the deviations of one 
of the variables and using this measurement to 
alter one of the parameters in such a way that 
the steady state values are restored. The flow 
rate of the coolant immediately suggests itself 
as the parameter to control but at this stage 
the question is left open whether the deviations 
of € or » should be used for this. There are three 
distinct types of perfect control to consider, in 
which the parameter deviation is made propor- 
tional to the value, the derivative or the integral 
of the variable deviation. Only perfect control, 
in which this is attained without time lags, will 
be considered in this paragraph, the consideration 
of imperfect control being the concern of par. 6. 
Moreover the possibility of control from the local 
linearization only will be studied, for if it is not 
possible to render a steady state stable locally 
it will certainly not be possible in the large. 
Let &,, n, be the concentration at a steady state 

and the parameters «, 8, y be defined as before 
and define the parameter 

aU ashe 

(: == rate of change of heat rejection 

\de/s with respect to variation of 
coolant flow rate. 


Using an asterisk to denote the deviations of any 
quantity from its steady state value the local 
linearization of equations (2.5) and (2.6) is 


(4.1) 


— 6q,* (4.2) 


g* = — (1 + a) &* + y7* 


yt = ag - ya" 

The method of control is to put ¢,* proportional 
to &* or »* or to one of their derivatives or 
integrals. The possibilities are as follows; in 
which for later simplicity the proportionality 
constant has been written K/5. 


Mode of 


Integral 
control 


Proportional Derivative 





(K/3) »* dt 


tn* 
(Kk §)° 


q-* (kK 5) 7* or £* 
dt 


q-* (K/5) 9* or ¢* (K/3) p »* (K/5) 9*/p 





In the second line of this table is given the 
relations between the Laplace transforms of the 
which are bar. The 


variables, denoted by a 


Laplace transform is defined by 
- 


e~?" &* (r) dr 


(4.3) 


and applying the Laplace transformation to the 
equations (4.1) and (4.2) there results 


(p 1 


x) f e. y "* 


/ 


(¢*) 


; * 
v) 9 


(4.4) 


(p 1 B (»*) 69,* (4.5) 


where (£*)° and (»*)° are the initial perturbations. 
If g, is uncontrolled, i.e. 7,.* = 0 the solution is 
A (p)é * (p 1 


A (p) 9* (p 1 


B — y) (*) 


‘ 


y (n*)° (4.6) 


x) (7*) x (£*) 


where 


A (p) p* d, p- d; (4.7) 


the determinant of the left hand side is a quadratic 
with (3.10). The 
relations gives the solution (3.9) and the condition 
for stability is again that the roots of (4.7) 
When gq, is 
controlled, one of the relations from the above 
table is inserted into (4.5) and the polynomial 


identical inversion of these 


should have negative real parts. 


A(p) is correspondingly modified. This will be 
done for the various modes of control mentioned 
above. 


(i) Proportional control on &. 


(K/8) &* the 


Here q.* 


and determinant A becomes 


p+-1+- al 4 
K +a p+1+s—y 


The for control is d, > 0, 
dy + yK > 0, so that if d, < 0 (i.e. the operating 
point S is in D, and D,) control can never be 
achieved. If d, > 0 however, it is only necessary 


A(p) = 


=p" +d,p +d, . yk 


condition stable 
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to make K positive and K > —d,/y for the 
controlled point to be stable. 

(ii) Proportional control on »*. 
(K/3) 9* and 


Here 7,* = 


A= p* + (d, + K)p + (d,; + K(1 + a)). 
Thus by making K > max { — d,, — d,/(1 + x)} 


both coeflicients can be made positive and stability 


eo” 


is assured, 

A similar analysis of derivative and integral 
modes shows that the following situation obtains 
with regard to local accessibility to control. 
A denotes that if the operating point S (a, 8, y) is 
in the region, control can be attained, N that it 
cannot. 


Mode Proportional Derivative Integral 





Region 


D, D, Dy D, Dy Ds D, Dy Ds 


NAN 


Control on &* 


Control on »* 





It is clear from this table that it is preferable 
to control on measurements of temperature 
deviation and that integral control is invariably 
too sluggish to succeed. 

Before passing on to imperfect control it is 
interesting to consider whether if all three modes 
are present the first two can be made to overcome 
the the third. For this, put 
q.* = (RK, K,,p) (9*/5) and obtain the 


cubic, 


slow ness of 


Ky, p 


A=(1+K,)p*+{d,+K, 4 
+{d,+ Ks 


K, (1 + %)} py 
K, (14 a)}p + K,(1 + 2). 


This has stable roots if K, > 0, (1 + K,) > 0, 
d,+ K, + K,(1 +) >0 and {d, + K, + K, 
(1 + «)} {dp + Ky + K,(1 + a)} > Ky (1 + Ky) 
(1 + a). This can be achieved by making K, and 
K, sufficiently large. On the other hand, if K, 
and Ky are fixed the first and last conditions put 
a bound on Ky, the amount of integral control 
that can be tolerated. 


5. IMPERFECT CONTROL 


A servo-mechanism so perfect that signals of 


temperature deviation produce proportional 
changes in the flow rate of the coolant without 
either delay or error does not exist. In practice 
a thermocouple reading will actuate a motor 
which will turn a valve and at all stages there will 
be errors and delays. Following the accepted 
practice of servo-mechanism theory, all these 
factors may be lumped into one transfer function, 
writing 

$9." 


: K G(p) 7*. (5.1) 


The way in which this transfer function is 
constructed is described in the texts (see, for 
example, CuestNut and Mayer [4]) and it is 
not our concern to discuss this here. It is impor- 
tant, though, to find out whether the time lag of 
an imperfect system which is ultimately a 
proportional controller can destroy the possibility 
of control. It is shown in CuestNut and Mayer 
that a transfer function ultimately giving pro- 
portional response can be of the form 

G(p) = M(p)/N (p) 


where 


M (p) = p”™ 
(p + #4) (p 


m 
m,p 


and 


N (p) = p" 
(p - Vv - v eee | v,) 


The », and », are reciprocals of certain time 
constants of the mechanism and are necessarily 
positive or the servo-mechanism itself would be 
unstable. Thus the coefficients m; and n,; satisfy 
the Routh-Hurwitz conditions 


m, 1 
m, > 0, 


| m1 0 | 
> 0, 


Ms My | Ms 


Ms My M, - 0, ete. (5.4) 


Mm, Ms 


Such a servo-mechanism inserted in our system 
would give stable the 
A(p) of the rational function 


control if numerator 


(pP+1+a){p+1+B—y7+KG(p)}+ay 
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has all its roots with negative real part. Now 


A(p) (p* + d.p 4 d,) N(p) 
+ K(p +1 + «) M(p) 


(d, + n,) p" 
d,) p” 
K{p™’' + [m, +1 
+ [m, + (1 x) m,] p” 
.(1 + a) m,,} 


eer 


+ (mg + d.n, dyn, 


x) p™ 


p' |, p' 


max (m+1,n-+4 2) If m>astl 

coetlicients will be of the form 
l= N+ KM, where M, = m,; + m;_, (1 + 2). 
From the Routh-Hurwitz determinant of order q 
a term A‘ « (the R.H. determinant of the M,) 
may be taken out and all other terms will be 


of order K*~' or less. For example, if m n + I, 


where | = 
all the 


N,1 
N,N, 


But the determinant of the M, is the R.H. 
the 
M(p) and so is positive. Hence by choosing a 
the R.H. determinants of 
the J, can be made positive and stability assured. 


determinant of polynomial (p+ 1+ 2) 
sufliciently large K, 


If m n. 


The factor of K*"' 1)" order R.H. 


determinant is 


in the (g 


1 0 0 
+ x)My 1 
(la)m,. 


m,+ (1 
m,-+ (1l+a)m, m,- 


+ajm, m,+ (1 
(l+a)my m+ 


which, since the first line is zero except for its 
first term, is simply the gth order R.H. determinant 
of (p + 1 + a) M(p) and so is positive. Again, a 
sufficiently large value of K will succeed in control. 


The conditions, m =n + 1 and m = n, may be 
interpreted to mean that the ultimately pro- 
portional control is initially derivative or pro- 
portional. It is, however, more likely that the 
response would initially be more sluggish, for 
example, corresponding to integral control, i.e., 
m =n — 1. In this case there will be limitations 
on the coeflicients n,; which involves the reciprocals 
of the time lags of the system. Since the first 
term in which A can appear to be of service in 
forcing it to be positive is the second, there is a 
first which must be 
satisfied by making the time lags »; 
that the sum of their reciprocals n, is greater 
than — d,. 


e 


condition, d n, > 0, 


€ 


' so small 


The second condition is 


K (n, + d l—a«a 


‘ m,) 
d,n,(n, + d,) 


(ny, My Ns) 


d, d, ~O (5.6) 
since K >O 
condition on the v; 


imposes a yet stronger 


However this condition 
1 


which 
i 


can always be satisfied by making the »; 
sufliciently small, for the term (nm, mz — ng) is of 
order — 3 in them and being a R.H. determinant 
is positive and can be made to dominate the 


is clear that A and the n, are 


expression. It , 


related in a much more complex fashion and 
control cannot be attained simply by increasing 
the gain AK. 
time lags make the coeflicient of A positive, e.g. 
d, if increasing K 


It is first necessary to insure that the 


in (5.6), ny > 1+ a+ m, 
is to be of any service. 

2 a much more serious 
the first 


In the case m n 


restriction appears for while again 


condition is d, + n, > 0, the second is 


K < (n, my — ng) + dn, (nm, + 4.) + d, dy 


which is actually a bound on the possible gain K. 
Thus for m < n — 2 stability can only be assured 
by making the time lags of the system sufficiently 
small and an increase in the gain K will actually 
induce instability. Two examples may be used to 
illustrate this. 


Let ( Y 
» + »>+ we 

G(p) = } Hy } Me 
(p + 1) (p + vg) 

be the transfer function of a servo-mechanism 

both initially and ultimately proportional. The 


characteristic polynomial 
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p* (vy + ¥ 
T p { Y, V2 
+ dy + K (pm, 
Pp [», vy d, + (vy V2) dy 
+K {i Hy + (1 + %) (my + 2) }] 
1% Me d, + K (1 + %) py Me , 
The first and fourth R.H. determinants give 


vy, + vw +d, + K >0 and 


y%4d,+K(1 x) fy fy > O 


which can be satisfied for sufficiently large A. 
The the 
conditions are 


leading terms in second and third 
K* (4, bog + 1 


k® (4, 


 «) and 


' 


fg+ 1 4 x) {os Mg + (1 %) (py Ha) 5 


so that by making K sufficiently large these 
conditions too can be fulfilled, however large the 
time lags may be (i.e., however small the v; and 
H;)- It that the 


coeflicients of K contain the sum and product 


is to be noticed, however, 


of yz, and py. Thus it will be much easier to attain 
control if the «, are large, i.c., the corresponding 
time lags small. 


In the case of an initially integral controller 


' PP 
(p + ¥,)(p 


Lp 


G(p) = (5.8) 


Al P) 


where 1, 


lp” + Isp 


Vy) d, + dy i 
Vy) dy K (1 + + ye) 
Ku (1 — x) 


VM, d, + (ry 





VM dy 
and the Routh—Hurwitz conditions are 


(5.10) 


(5.11) 


(5.12) 


Il, = rd, + Ku(l +2) >0 (5.13) 


For any given set of values A, yu, «, 8, y these 


inequalities define a region of the »,, v, plane 
within which the point N(»,, v,) characterizing the 
time lags of the system must lie. Fig. 4 shows this : 
(5.10) is a straight line bounding N away from 
the origin, (5.13) is a rectangular hyperbola and 
(5.11) and (5.12) give curves of similar shape to 
the hyperbola of which only one is shown. The 
region is symmetrical about the line », 
it must include at least the region defined by 


v, and 


-d,+1+a+4, 


¥) ve > Ky (1 (5.14) 


x) (— dy) 








Ky(l+a) de 


-~de+1+a+ 
-% de atu 


4. Restrictions on the time lags of an imperfect 
proportional controller. 


6 A MeraAsuRE OF LocaL CONTROL 


It is convenient to introduce a numerical measure for the 
degree of stability of a steady state and hence of the 
effectiveness of control. Suppose that W (f) is the response 
of one of the system variables to a unit disturbance in 
one of the others, then 

fe) 


W(t)? dt 
* v0 


(6.1) 


is a measure of the total disturbance caused by this 
impulse, For example if in the uncontrolled equation one 
puts (»*) 

€ and » are 


1 the responses We» and W,,. of the variables 


Weft) | : ) (eon #2") 
w) Wwe 


Walt) 


Y 9 
| | (cy + 2 + adel’ 
bod | a9 


If and only if w,; and wg have negative real part the 
integrals J will converge and 


(we + x) 2") 
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Tg y 2 d, dy 


Ton id + (1 + a)?} 2d,d, 


The stability of the steady state is thus greater or less 
according as the product d, dy is large or small and this 
product gives a useful measure. Further, the attainment 
of control I’ may be defined by 


(d, dy) uncontrolled 


r= 1 (6.2) 


(d, dy) controlled 


a quantity which is zero for no control and unity for 
perfect control. For example, for the control of ¢,* on »* 


I’ has the physical meaning that it is the difference 
between the total disturbance for controlled and un- 
controlled states as a fraction of the uncontrolled distur- 
bance. For an originally unstable state this meaning is 
lost since the total disturbance is then infinite but the 
measure defined above in (6.2) retains a meaning for 
sufficiently good control approaching unity as the control 
is increased. A measure of control on similar lines to this 
has been introduced by Bitous, Brock and Pirer [3] 
in their treatment of the locally linearized case. 


CONCLUSIONS 


In this paper the application of the method of 
local linearization to the analysis of chemical 
reactor control has been discussed. It may be 
used to determine the stability of the immediate 
neighbourhood of a steady state, and hence the 
possibility of control by various means. For 
clearly, if stability is not attained in such a 
neighbourhood, there is no hope of control in any 
larger sense. It is found preferable to measure 
the deviation of the reactor temperature from 
that of the steady state and use this signal to 
control the rate of cooling. If the response of a 
perfect controller is proportional to this deviation 
or to its derivative, control may always be 
attained by sufficient amplification of the signal. 
A response proportional to the integral of the 
deviation, however, is too sluggish to be effective. 
With imperfect control, when there are certain 
time lags in the system, it has been shown that, 
even when the response is ultimately proportional 
to the deviation, these time lags must be kept 
below certain bounds. 
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The evolution of proportional control 
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Abstract 
reactor temperature and uses them to give a perfect proportional control of the rate of cooling. 


In this part it is assumed that the control system receives signals of the deviation of 


It has already been shown that control cannot be attained until the proportionality constant 
is great enough but it appears here that even when this minimum value has been exceeded, 
the control may only extend to comparatively small disturbances. To see how this comes about, 
it is necessary to follow the behaviour of the reactor through the variation of the control parameter, 
and it is shown that the methods of non-linear mechanics make this possible without excessive 
labour. There comes a point, however, when some computation is necessary and this is ex- 
emplified in a complete description of a simple system, 


Résumé— Dans cet article les auteurs supposent que lappareil de controle recoit les signaux des 
variations def du réacteur et les utilise pour donner un contréle proportionnel, parfait de la vitesse 
de refroidissement. On a déja montre que le controle ne peut étre obtenu tant que la constante de 
proportionnalité n'est pas suflisamment élevée. Mais il semble ici que méme lorsque la valeur 
critique n'a pas été dépassée, le contréle ne peut s‘étendre qu’a des perturbations relativement 
faibles. 


variation du paramétre de controle ; 


Pour verifier ce comportement, il est nécessaire de suivre la marche du réacteur pendant la 
les auteurs ont montré que les méthodes & mécanisme non 
linéaire, rendent ceci possible sans trop de travail. Mais il peut arriver que certains calculs soient 
nécessaires, et ceci est illustré dans une description compléte du systéme simple. 


Zusammenfassung—In dieser Arbeit ist angenommen, dass das Regelsystem Signale iiber die 
Abweichung der Reaktortemperatur empfingt und hiernach eine vollstindige Proportional- 
regelung der Kishlgeschwindigkeit durchfihrt. Es ist schon gezeigt worden, dass eine Regelung 
nicht kann, ohne dass die Proportionalitdtskonstante gross genug ist; hier 


scheint sogar nach Uberschreiten dieses Minimalwertes eine Regelung nur bis zu vergleichsweise 


erzielt werden 


kleinen Stérungen ausfihrbar zu sein. Um zu sehen, wie es dazu kommt, muss man das 
Verhalten des Reaktors bei Variation des Regelungsparameters verfolgen und es zeigt sich, dass 
die Methoden der nichtlinearen Mechanik dieses ohne allzu grossen Rechenaufwand gestatten. 
Indessen geht es nicht ohne Rechnung ab und diese wird an der vollstindigen Beschreibung 
eines einfachen Systems gezeigt. 


Heat 


discussion so 


7. Tre TRANSFER FUNCTION 
Tue far itself 


side the reactor and », and q,, the inlet tempera- 


has concerned ture and flow rate of the coolant, and this might 


wholly with the local properties of the functions 
involved and before going on to more general 
considerations in the large, the heat transfer 
function demands a little consideration. In fact, 
all that is required is a function expressing the rate 
of heat rejection in terms of » the temperature in- 


well be an empirically determined one. Without 
implying any restriction on the method here pre- 
sented it will be convenient to derive a fairly simple 
approximate expression for U. In the first place, 
it will be supposed that it is proportional to the 
difference between the reactor temperature and 


*For ease of reference, the numbering of the sections and figures is continuous with Part I. 


**Present address : 


The Mathematical Institute, University of Edinburgh. 
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and the mean coolant temperature, 9, and, 


following respected correlations, to 
Thus 


secondly, 
the four-fifth power of the flow rate. 


U(n) = A (4) Me) (7.1) 


Again, the total rise in temperature, is 2 (%, ".)s 
so that since the heat rejected goes to warm up 
the coolant 


U(») 2B (7) (a. ”,) (7. 


Eliminating 4, from these two equations 
O (93 ne Ge) 
q. | q. 7 -~ «€ 
AB (7 7 A+B 7.3 
| ‘) 1— Ne) | ° (7.3) 


A is a constant involving the cooling area and 
materials of the heat exchanger, B involves the 
heat capacity of the coolant. 

For the perfect proportional control assumed 
here, the deviation of g, from its steady state 
value g,, is to be made proportional to the 
deviation of temperature » from the desired 
steady state value »,. Naturally this must be 
within the 
and it 


range of movement of a_ control 
is assumed for simplicity that 


half its 


valve 
in the steady state the flow rate is 


maximum. Thus 


4% - Ves 
q 


K (» (7.4) 


»,) for I» 


where » 


r 


M (‘*) is the half range of tempera- 
For 


ture which causes the valve to move. 


while for 
(7.6) 


Introducing (7.4) into the expression (7.3) for 
U and expanding in a Taylor series 


ABg,, 
U(») F lt 


+ B(q., 53 "e) 
+ B(q., 9)" 
$ B(q., 9)" 


KAB | 4 
(4 





(y — 4.) (9 — ,) 


when In, _— | < 7, 
now be written with convenient simplicity, will 


This expression, which can 


be of good accuracy since the (g,/¢)°* in the 
denominator is of comparatively little effect. 
Therefore 


U,,A(y : Ne)s 
Udn—7-) + k(n 
0, 


U(») 
f (7.8) 


k »,, the function U depends 


when 





Since U,, — U 
on the proportionality constant & in a peculiar 


way, which is perhaps best illustrated by Fig. 5. 
For very high gains, k, this reduces to off-on 
control in which the slightest drop in tempera- 
ture shuts off the coolant and a slight rise opens 
the valve fully. This represents the limit of this 
form of control and will be considered in par. 13. 
In the calculations which follow it will be assumed 
that the range of the valve is such as to make 
the quadratic expression (7.8) valid. 











Fic. 5. The heat transfer function. 


8. Tue Puase PLANE 

The notion of the phase plane has long been 
in use in non-linear mechanics for equations of the 
type (3.1). It allows one to get a picture of the 
behaviour of the system in the large, in contrast 
to the local linearization which gives information 
only in the neighbourhood of the critical point. 
Since time does not occur explicitly in equations 
(3.1), » may be taken as the independent variable 
so that 
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& _ Pl, 0) 
” Q(é, ”) 


(8.1) 


which is a single non-linear equation for € as a 
function of ». The value of dé/dy is unique at 
every point where P and @ do not simultaneously 
vanish ; such points are called ordinary points. 
The values of £, » where P=Q=0 are, of 
course the steady states already discussed and 
these are known as critical points of the phase 
plane. Starting from an ordinary point the 
equation (8.1) can be integrated and gives a 
curve in the phase plane, a trajectory, which 
represents the life history of the system. Inte- 
grating the pair of equations (3.1) and plotting 
& versus 7 would give the same curve but now a 
time scale could be marked on the trajectory. 
This may be illustrated from the series of figures 
derived for the special case discussed later. For 
example, consider the trajectory DA in Fig. 11 
starting from the point £ = 0, » = 1-75. This 
represents the history of a reactor when initially 
the temperature is 1-75 and concentration of X 
zero and it is seen that the temperature rises 
only slightly while the concentration of XY goes 
up to some 96 per cent of its inlet value and at 
this critical point the steady state is reached. To 
see how quickly this comes about one can compare 
the curves marked A in Fig. 10 (a) and (b) where 
the concentration and temperature are separately 
plotted against time. It is seen that the approach 
is quite rapid and smooth as is characteristic of a 
nodal critical point. 

The phase plane representation enables one to 
make considerable qualitative predictions about 
the behaviour of the system without integrating 
the equations in detail. Much 
so-called topological method will be made in par. 9 
but it will be useful here to elucidate certain 
general points. It is first noticed that a trajectory 
must always remain within a finite part of the 
phase plane. For low temperatures P +0 and 
U <0 so that by (2.6) 7 >0 while for high 
temperatures P and U are large and 7 < 0. 
Thus at low temperatures the temperature always 
tends to increase and at high temperatures to 
decrease. Moreover when £ = 1, £ < 0 and when 
é=0, P=0, and é=1>0 so that from the 


use of these 


upper and lower ¢ boundaries a trajectory always 
tends to move inwards. By the second theorem 
of Bendixson (see Minorsky [10]* p. 78) it can 
be stated that any trajectory 


either (i) 
or (ii) 
or (iil) 


approaches a critical point 
approaches a closed curve 
is itself a closed curve. 


An example of (i) occurs in Fig. 11; Fig. 17, 
corresponding to a different value of the control 
parameter, illustrates (ii) and (iii), Here the 
trajectories starting from the points A and C 
wind around and gradually approach the closed 
curve DEF (shown as a broken line) the one from 
and the other from the The 
closed curve DEF is itself a trajectory and is 
called a limit cycle. It represents a sustained 
oscillation of the temperature and composition 
and under the conditions represented by this 
figure the system will always tend to get into 
this state. 


inside outside. 


In a previous paragraph it was seen that 
through a saddle point there are four trajectories, 
two which approach it and two which diverge 
from it. The two which approach it must come 
from somewhere, either from an unstable critical 
point or limit cycle or from the boundary. For 
example, in Fig. 11, it is found by integrating 
backwards that these trajectories come from 
points F and G on the lower and upper boundaries. 
On the other hand the trajectories going outwards 
must take one of the three alternatives given 
above and in Fig. 11 one goes to the node at A 
and the other to the focus at C. As was pointed 
out previously [2], the separatrix /' BG is of great 
importance for starting up the system. It divides 
the phase plane into two parts such that from 
any initial conditions represented by a point in 
the left half the low temperature steady state is 
approached, while from any conditions in the 
right half the high temperature steady state is 
attained. 


9% Tue Puase PLaANne witu ContTROL 
It is interesting to see what can be said about 
control without doing any detailed calculations. 


*All references are listed at the end of Part I of this 
paper (p. 131). 
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It will be seen that a very good general picture can 
be built up from quite general arguments about 
the disposition of the critical points in the phase 
plane. Such arguments are well known in non- 
linear mechanics under the name of topological 
arguments. This study, which is further aug- 
mented by some simple calculations in the next 
paragraph, is both valuable in itself and a 
necessary preliminary to the more detailed calcula- 
tions described later. 

If the expression for U derived in the previous 
paragraph is inserted into equation (2.6), the 
pair of equations for £ and » are 


€=1—£€— P(é, n) 


n = (1 + U,) (9 — 9) — k(n — 4.) (9 — 1) 


where 9 = (mn) + U,7,)/(1 + U,). 
previous procedure £ may be eliminated between 
the two equations giving 


(1 + U,)(n — 9) + My — 9.) (9 — 4) = 
P (£,(»), ) 


where the right hand side is again the curve 
R in Fig. 1. With k = 0 i.e. in the absence of 
control, the left hand side is the straight line L 
but with any k > 0 it is a quadratic curve 
intersecting L at » = », and » = »,, and increas- 
ing in curvature with increasing k. It is the 
intersection of these new curves, e.g. Q and Q’, 
with R which determines the new steady states 
for the various k. 

In the first place it is clear that unless », is a 
steady state temperature for the uncontrolled 
system it will not be the temperature at which 
the reactor will be controlled for all values of 
k, for the increasing curvature of Q will cause it 
to intersect R at a varying point always different 
from », unless L and R intersect there. If, 
however, », is the temperature of a steady state 
B of the uncontrolled system it will always be 
one for all k. The other critical points will move 
as for example when Q’ cuts R at A’ and C’ but 
ultimately a sufficiently large k may give a 
single intersection such as is shown by the curve 
Q. Then there is some hope of having attained 
control in the large sense as well as locally. 


(9.1) 


(9.2) 


Following the 


(9.3) 


To discover the position and nature of the 
critical points for all k it is convenient to reverse 
the procedure and determine the value of k 
which gives a critical point at a chosen tempera- 
ture. From é = 0, the relation € = €,(y), which 
must be satisfied for a critical point, may be 
obtained, and hence from (9.3) the value of k 
for which a critical point is at &(n), 7, namely 


— P(é, a) + (1 + U,) (4 — 1) 


k= me 9.4) 
(7 — %)(m — Ms) { 


Even if the more accurate expression (7.3) were 
used for U it would still be possible to solve for k 
without much trouble. Now that € and » are 
known so are the parameters a, 8, y and d,, d, 
associated with the critical point and so its 
nature and stability are determined. 

To illustrate the procedure a simple example 
is introduced, which despite its apparent simplicity 
shows the quite intricate evolution of reactor 
control, Although the simple character of the 
example allows performance of some of the steps 
more easily, the method outlined here is not in 
any way restricted by the complexity of the rate 
equation, P(£,), and its adaptation will be 
immediately obvious to the reader. The example 
concerns a first order irreversible exothermic 
reaction, € being the concentration of the reactant 
X being destroyed. Since (— AH) was defined 
as the heat of reaction for the formation of X it 
must be taken as negative in the second equation 
thus changing the sign of the reaction term. Thus 


P(E, n) = € 


exp 30 (3 a a} =£€¥Y (9.5) 


J 


; 
2 n) 


and the equations are 
§=1—¢-—€Y = P(£, ») 


7 = — An — 9) — Ky — 9) (9 — 2) + EY 
= QE, ) (9.7) 


If, = no = } = 1-75 and yn, = 2 these equations 
have a critical point at §=4, » =2. To fix 
ideas still further this may be regarded as a 
reaction taking place in a _ vessel of 501. 
capacity. If the flow rate g is 50 1/min the hold- 
up time or unit of 7 time is 1 min. If the activa- 
tion energy of the reaction is 20,000 cal/mole, 


(9.6) 
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the unit of » corresponds to 200°K. The critical 
point of interest is thus at a temperature of 400°K 
or 127°C and with the reactant concentration in 
the tank at 50 per cent of its inlet value. Most 
of the concern is with the temperature range 


vd ~~ 


n = 1-75-2-25, i.e., from 27 to 227°C. 


The parameters of the steady states are 


Y; 8 1 k(2y - ”,)3 


50 £ Y /y* (9.8) 


y 
At the point » = 2, where an attempt is being 
made to get control of the originally unstable 


steady state, x = 1, 8 = 1+ k 4, y = 625 and 


d, = (k-9)/4 dy (2k-9) 4 (9.9) 


Hence there is no hope of control until k > 9 
when both d, and d, become positive and the 


point stable. Solving the quadratic, w* +d, 


0, at this point there are two roots 


b(k—9) + £y(k 


Thus the character of the controlled point is as 


| dy 


5) (k—45) (9.10) 


1.2 


given in the following table. 





Range of k Nature Stability 
Unstable 
Unstable 
Unstable 
Stable 


Stable 


Saddle-point 
Node 

Focus 

Focus 

Node 





The technique just outlined for determining 
the moving critical points may be applied for by 
(9-6-9-8) 


E(m)=1/(14+ VY) (9.11) 








| 
| 
3D 325 


‘75 180 ~—SOST SD 20 - 


Temperature, 7 


35— 


Fic. 6, Position of critical points for varying k. 


' (9.13) 





4 


From these the position of the critical points 
is obtained and by solving the quadratic for w 
the characteristic roots their 
determined. The table summarizes the 
information obtained from this calculation. 

In Fig. 6 is shown the position of these critical 


and nature is 


below 


points in the phase plane, the value of & being 
inscribed by each marked point and the letter 
corresponding to the above table : all, of course, 
fall on the curve € = €,(y). In Fig. 7 the value 
of w, and w,, the two characteristic roots, is 
shown and the letters on these curves show the 
parts of them which correspond to the stages in 





Low tem perature 
Range of k Nature Stability 
Stable 
Stable 
Stable 
Stable 
Stable 
Stable 
Stable 


0 Node 
O-O-34 Node 
O34 —-O-847 Focus 
OS47 3-33 Focus 
3-33 -3-64 Node 
364 -4:55 “F Node 
45 5-00 Node 
5-00 


High tem perature 
Stability 


Critical point 
Nature 


Stable 
Stable 
Stable 
Unstable 
Unstable 
Unstable 
Unstable 


Focus 

Focus 

Focus 

Focus 

Focus 

Node 
Saddle-point 

Coalesence 
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the above table. We may now inquire what this 
knowledge together with general topological 
considerations allows us to say about the phase 


plane as a whole. 





Characteristic root 


, Aon temp. critical point 
Real part 
Imaginary part 
Mean of real roots — — — — 





a 
k— proportionality constant 


Fic. 7. Characteristic roots for varying k. 


From earlier work there is some knowledge of 
the phase plane when k = 0. From the saddle- 
point there issue two trajectories which go to 
the stable critical points and from the boundary 
there arrive two trajectories which divide the 
phase plane into two parts, all trajectories in 
the one part going to the stable critical point 
there situated. When k increases beyond 0-34 and 
the low temperature point becomes a focus no 
significant change occurs for the trajectory from 
the saddle-point still proceeds there. But when 
k increases beyond 0-847 a radically different 
situation exists for then the high temperature 
singularity has become unstable and the trajec- 
tory leaving the saddle-point can certainly not 
arrive there. Three topological possibilities exist 
which would take care of this situation; these 
are shown in Fig. 8. Either the trajectory 
arriving at the saddle-point, B, now comes from 
the unstable critical point, C, and the trajectory 
leaving B, which formerly went to C, takes a 
detour around the bottom and goes to A, as is 
shown in Fig. 8 (i); or the two arms of the 
separatrix join together to form a closed loop 
and the trajectories within this loop tend out- 


wards toward it while those without skirt it and 
go to A, as is shown in Fig. 8 (ii). Or a stable 
limit cycle, L, grows around C and attracts both 
the trajectory issuing from B and also all trajec- 
tories to the right of the separatrix, FBG. In 
the last case the division of the phase plane into 
halves is maintained, though the stable critical 
point of the right half is replaced by a stable 
limit cycle. In the second case all trajectories 
starting outside the closed loop will eventually 
arrive at A the low temperature state, while 
the loop itself forms a pseudo-limit cycle to 
which the trajectories from within tend, but 


since B is a critical point where £ and » vanish 


the system will spend increasingly long periods 
in the neighbourhood of this point but never 
ceases to make excursions round the loop. In 
the first case, however, all trajectories, except 
the separatrices, tend to the one stable critical 
point A. It will later appear that the first 
alternative is what occurs but it is not possible 
to decide between them at this stage. 


























F\ 


Fic. 8. Topological possibilities when k > 0-847. 
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We do, however, know that as & increases 
through 4-5 the critical point C passes through 
the fixed point B and after this it is the moving 
critical point which is the saddle-point and the 
fixed point B which is an unstable node or focus. 
Moreover, one can determine the directions of the 
separatrices through the fixed saddle-point B, 
following the reasoning of par. 3, and they can 
be shown to vary only slightly. At & = 4-5 the 
separatrices pass through £ = 4, » = 2 in the 
directions d&/dyn = — 0-5 and dé/d» 0-32. 
If the third were the correct alternative of those 
presented in Fig. 8 the limit cycle would be forced 
to assume a very pointed vertex near k = 45 
and so a break-down into the form of (i) or (ii) 
would seem inevitable. But the alternative (ii) 
is unsatisfactory at this point also for a closed 
curve cannot be a trajectory without containing 
a critical point and if C is made to coincide with 
B it would leave the loop in (ii) without any 
critical point within it. The same objection cannot 
be raised against the configuration of Fig. 8 (i), 
however, for the closed loop when C coincides 
with B would here be not one trajectory, but two. 
Thus it already appears that (i) is the correct 
alternative, though the detailed calculation for 
confirmation will be determining. As *& increases 
from 4-5 to 5-90 it is the saddle point which is 
now moving with the fixed unstable point within 
it, but this is basically the same configuration as 
before. At k = 5-90, however, the node A and 
the saddle-point C coincide and that part of the 
loop going in anti-clockwise direction vanishes 
while the trajectory leaving the fixed unstable 
point B and going to the saddle-point C can no 
longer arrive there. Thus a limit cycle is left 
behind which is traversed in a clockwise sense and 
all trajectories from the still unstable point B wind 
outwards towards it. For values of k between 
5-90 and 9, the reactor cannot but behave in a 
highly unsatisfactory manner, always ultimately 
tending to a regular cycle of temperature and 
concentration variation. Without full calculation 
nothing further can be said about the behaviour 
of this limit cycle save that ask -> 9 it can either 
close down on the fixed point = }, » = 2 and 
vanish as this becomes stable or it can persist 
after k = 9. In the latter case an unstable limit 


cycle must grow up within it surrounding the 
now stable critical point. The method of the 
following paragraph enables us to decide which 
alternative is true. 


ATk=9 


The study of the interesting phenomena which 
can occur when a critical point changes its 
stability was initiated by Porncare [11] and is 
known generally in the texts on non-linear 
mechanics as “bifurcation theory”. It goes 
essentially beyond what can be studied by the 
linearized equations though it is limited in that 
it can only say what takes place in the immediate 
neighbourhood of the change. It will, however, 
allow us to say which of the alternatives stated 
at the end of the previous paragraph is the correct 
The main ideas of the method will be 
first outlined, and then the details of the method 
given for the general case, and finally, the example 
being studied will be worked out. 

Transition from instability to stability will 
normally take place when the real part of the 
characteristic roots w of a focus turn from 
positive to negative. At the transition the 
roots are wholly imaginary and the solutions of 
the linearized equations give a system of concen- 
tric ellipses in the phase plane about the critical 
point. By a suitable linear transformation of 
co-ordinates these ellipses can be turned into 
circles, and the equations into an equation for r, 
the distance from the critical point in the trans- 
formed phase plane as a function of @, the polar 
angle and the parameter k. This equation takes 
the form 


dr 
d@ 


10. Tue TRANSITION 


one. 


= R,(0,k)r + R,(0,k)r? +... (10.1) 
and for the transition value of k (say k*), 
R, = 0. This is to be expected, for the linearized 
equation, i.e., that obtained by neglecting R,, etc., 
is to give a system of circles for which dr/d@ = 0. 
If ry is the value of r when @ = 0, an attempt is 
made to find a solution of the form 


‘= r(O,k) = rou, (O,k) + rot, (O,k) T see 


(10.2) 


by substituting and comparing the powers of r,. 
This will represent a closed trajectory if after 
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the angle @ has increased by 27 the same value 
of r is obtained, i.e., if 


r—ro=Tof Uy (27,k)—1} +17 uyg(2a,k) +... (10.3) 


is zero, Again at the transition value k* the 
coefficient of r, must vanish and it is also found 
that u,(27,k*) = 0. Expanding (10-3) as a 
Taylor series in powers of (k — k*) and retaining 
only the lowest order terms the condition for a 
limit cycle is 


ik — ae) 


aad) 4+ r2ug(2m,k*) =0 (10.4) 
dk knk® 


Thus when the first term is of the opposite sign 
to u, there are real values of r, different from zero 
which give limit cycles. 

At a transition from instability to stability 
[du,/dk}, »e< 0 hence if us >0 limit cycles 
exist for k > k* and these must be unstable since 
they surround a stable critical point. If u,< 0 
the limit cycles are stable surrounding the yet 
unstable critical point for k< k*, and they 
disappear for k > k*. It is in this way that the 
two alternatives may be distinguished. 


To formulate this programme in more detail and in a 
way which is immediately applicable to our equations, 
the following equations are considered. Let A = k — k* 
so that the transition of stability occurs at A = 0, then 
the roots of the characteristic equation (3.10) are 


where np (0) =0 (10.6) 


Let the equations be 


& = P(é,n) = a(A)(é — &,) + OC) (n — 1%) 


+ p*(£, 9) (10.7) 


QE. 9) = (A (E — &) + (A) (y — 29) 


+ q* (€, ) (10.8) 


and the local linearization about the steady state €,, 7% 
is obtained by omitting the functions p* and q*. The 
characteristic equation of which the roots are (10.5) is thus 


w*® — {a(A) + dA)} w + {a(A) d(A) — b(A) e(A)} = 0 (10.9) 


By making the substitution 
é— &, =2(u¢— a)u + 2vv (10.10) 
7—-% = —2cu (10.11) 


the canonical form of the equations (10.7) and (10.8) 
results 


“= pu — wo + p(u,v) (10.12) 


D = vu + po + q(u,v) (10.13) 


where — 2ep = 


2v cq 


— q* (u,v) and 


cp* (u,v) + (uw — a)q* (u,v) (10.14) 


These functions p and q contain the non-linear part of the 
equations and when expanded would contain only squares 
and higher powers of u and v. Though the formulae do 
not explicitly show it, it should be remembered that ,, », 
etc. are still functions of the parameter A. 


Now introduce polar co-ordinates in the u, v plane 


u=reos@, v=rsing (10.15) 


with the result that equations (10.12) and (10.13) become 


dr P 
r— = mw (A)r* + rcos@ p(r, 0) + rsin@q(r, 8) 


a (10.16) 


_ do : 
2— = Ayr? 4 


r ai rcos 0 qr, @) — rsin @ p(r, @) (10.17) 


and dividing one equation by the other 


dr ” pniArr + pcos @ +qsin@ 


= - —$— 10.18 
do v(A)r + qceos@ — psin#g ( 


this equation is expanded into the form 


R,(@, A)r + Rg (0, A) r® + RAO, A) PF +... (10.19) 


p (r, 0) = 1? (peg cos? 8 + py, cos @ sin 8 + pog sin® A) 
+ 1 (pao cos® 0 + pa, cos* @ sind +... 


and q is of similar form, then with rj; = 
Sij "Vij - # Pi: 


Ry(O, A) = wl A)/oA) 


vPij — HY and 
(10.20) 


1 
RO, A) = 3 [rao cos® @ 4 (711 + Seq) cos” @ sin @ 
is 


+ (Tog + 813) cos O sin® 8 + sg sin® @] (10.21) 


R;(0, A) = = lee cos* @ + (rg + Sq) cos® 8 sin 6 
v 
tT (ye + 81) cos @ sin® @ 


+ (Tog + 812) cos @ sin® @ + sog sin* @] (10.22) 


. [420 cos® @ 4 (411 — Peo) cos” @ sin @ 


+ (og — P11) C08 O sin? @ — pog sin® 0] Ra(O, A) 


If r =r», when @ = 0, a solution may be formed 


r(@, A) Tj 4,(0, A) + te" ug(@, A) 


, tee ug (0, A) +... (10.23) 
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and substituting and comparing coefficients of rg gives a 
series of recurrent differential equations 


du, 


i Ry(@, A) 


(10.24) 


dug 


- tg RO, A 
do fo M4 / 


u,* RAO, 4) (10.25) 


dug 


7 tg RYO, A) 


+ 2my Wg RelO, A) 4 uy* Rg, A) (10.26) 


By assumption r 
0 when @ — 0. 


exp | s e) 
' 


and the solution of an equation of the form 


etc. 


3% 


ro when @ 0 se 


that m4 


solution to (10.24) 


m,(@, A) (10. 


du Mm 
do ’ 


' Sie) 


d) - f(b) ded (10.28) 


it is necessary that 


| Ugl 2, 
| 

Ug(2r (10.20) 
linearized 
When 


0, for using 


When A 
equations for pf A) 


0 this is clearly realized for the 
0 and exp (27, /¥) I 0 
so that my 


(10.28) 


1 it is observed that uy (27, 0) 


Me (2.0) Reid) dé 


7 0 
and by (10.21) Ry is an homogeneous polynomial of odd 
degree in cos @ and sin @ so that its integral from 0 to 2 
vanishes. Again applying the formula (10.28) to the 
equation for ug and using (10.25) 


” 
2a 


2a P 
Ng(27, 0) [matz>. 0)| ‘ | Raid odd 
o 
« O 


and the first term on the right hand size is zero. From 


the well-known values of integrals of the form 
22 
cos*™ ¢ sin®” é6 dé the value of Hgli27, 0) can be culeu- 
rn 
lated from (10.22) in terms of the coefficients p,;; and ¢,;. 
If A # 0, various zeros of (10.29) may be expanded as 


Taylor series in A and this condition may be written as 


0 A fae 


—_— Mol2m, A) ir 
(0) . |, e a 


e 
lug(27, 0) ry” 


If this is to be satisfied the first and third terms must be 
0 (r4°) 
and the second term is negligible by comparison when A 


equal and opposite, ie. A = 0 (r9”) and so Arg 
is small. The size of the limit eyele for small A is thus 
given by 

phi) 


(10.30) 
WO) Ug(27, 0) 


If we (0) - 


on 
Phen rp 


0 the transition is from instability to stability. 
is real if either . 
(a) ug > 0, A> O; when an unstable limit cycle grows 


around the stable point, or 


(b) Wg <— OLA; 
unstable point shrinks around it and vanishes as it 


0; when a stable limit evcele around the 


becomes stable. 


For »’ (0) > 0 the transition is from stability to instability 
and either an unstable limit evele shrinks on the critical 
~ OLA 


point (uM, 0) or a stable limit evele grows around 


the unstable point (wg < 0, A > 0). 
Applying this technique to our example gives 
k — 9, A/8, 
VA 
2a pw’ (0) v (0) 


4) (36—A)/8 (10.31) 


7r/6<— 0 
u,(27, 0) = 78-0707 7 > 0 


and hence 0-0462 A" * forA >0O (10.32) 


Stabie limit cycle 


Stable critco! os 


The situation immediately after local control has 


Fic. 9. 
been attained. 


We thus arrive at the very important conclusion 
that the limit cycles do not vanish as & increases 
through 9 but rather that an unstable limit cycle 
around the controlled point within the remaining 
stable limit cycle. The situation in the phase plane 
for k slightly greater than 9 is thus as shown in 
Fig. 9. The circle of radius. r, transforms back 
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into an ellipse in the £, » plane within the nearly 
elliptical stable limit cycle. A trajectory starting 
from any point within the unstable limit cycle 
will spiral inwards to the controlled point but 
any other initial point will give a trajectory which 
tends to wrap itself around the stable limit cycle. 
Thus control has really only been achieved over 
the interior of the unstable limit cycle and a 
perturbation of the system which leads outside 
of this will result in the system performing the 
limit cycle. As & increases still further 
might expect the stable limit cycle to shrink and 
the unstable limit cycle to grow until they coalesce 
and evaporate leaving one stable critical point 
in the phase plane and control has finally been 
established. 


one 


This behaviour, however, cannot in 
general be fully decided without detailed calcula- 
It will be noted that all that 
required in the foregoing discussion is the expan- 
sion of the non-linear functions P( , ») and Q(€, ») 
as far as the third-order terms; this, though it 


tion. has been 


involves a certain amount of algebraic labour, does 
not demand the machinery required by the full 
scale integration of the equations and is, moreover, 
a valuable preliminary to the latter. 


In a certain case, however, a little more can be said 
about the behaviour of the limit cycle when k > 9 and 
since this applies to the present example the description 
of results, which follows in the next paragraph, will be 
0 the 
stable limit cycle is already of very small size and has been 
decreasing for increasing negative A, then it 
expected that the expanding unstable limit cycle will 


anticipated to explain this. Supposing that at A 
may be 


coincide with the stable one at some small value of A > 0 
and it should be possible to estimate this. 

If the stable limit cycle is small it will be nearly elliptical 
and its form will satisfy the linearized differential equations. 
If it is given by the equation 


¢2 eat 
Sé , «207 


then substituting from the linearized equations the co- 
efficients f, g, h must satisfy 


af + 0 
bof (a 


hypo t (10.33) 


djg+th-O 
bg + dh — 0 


and the condition for consistency is that the determinant 
d)(ad — be) 0. 
0 at a transition. 


(10.34) 


of these equations should vanish, i.e. (a 
This is satisfied since a 
The ellipse is thus 


d = 2y 


c &*7 4+ 2d én —bn* = d/g (10.35) 


The 
remembering that a + d 
found that it gives 


transformation (10.10), (10.11) is applied to this 


-0, v* = ad — be and it is 


‘ e ‘ 
24 yp? = d/4v* cg = 1,” 


“ 0 


and the expression on the right hand side is just the 
square of the radius of the circular limit cycle in the 


u, v plane, Le. 


(10.36) 


On the other hand, the reciprocals of the squares of the 
semi axes of the ellipse (10.35) are the roots of the quadratic 


” 
o- 
— 


d? 


b) 


1g [(c b) 


2 


b> — 47] 


Vile (10.37) 
If B is the semi-minor axis on eliminating g/d between 
(10.36) and (10.37) 


Be Py 8 ve [(c b) Vv (c b)* 4 | (10.38) 
which gives an expression for the growth of B when To is 
substituted from (10.32). In this case B = 0-052 A1/2, 
Anticipating a little the results to be exhibited in the 
next paragraph and referring to Fig. 21, it is seen that as 
k -» 9, i.e. A» 0, the minor axis of the limit cycle becomes 


quite small and is rapidly decreasing. In fact when k = 9 


A 0, BF 0-026, dBt/d\ = 0-180 


It is expected that the minor axis of the stable limit cycle 
for values of k near 9 will be approximately 


B* — 0-026 — 0-180 A 


0-052 A? both limit cycles will 
coincide and vanish, This gives a quadratic for A/? and 
0-025. 


and when this equals B 
from the positive root of this, A 


Thus for as low a value as k = 9-025 it is 
expected the limit cycles will disappear. It is 
perhaps regrettable that this state of partial 
control should have disappeared so soon in the 
particular example chosen, but there was no 
way of knowing this in advance. The importance 
of the phenomenon should not be overlooked 
for this reason, for in other cases it might be that 
the area over which control had really been 
obtained continued to be quite small long after 
the control parameter had passed through its 
critical value. Any imperfection of control, such 
as must certainly occur in practice, may be 
expected to the likelihood of this 
happening. 


increase 
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‘1G. 10. Concentration and temperature as functions of time. 


1l. DescriIpTion oF RESULTS 


An account of the results of extensive calcula- 
tions which have been done on the two equations 
(9.6) and (9.7) will now be given. This will give a 
complete picture of the evolution of control and 
show to what extent the predictions of the last 
two paragraphs have been valid. Those details 
of the calculation which are of general application 
will be discussed in the following paragraph ; 
suffice it to say here that the curves of Figs. 10-21 
are based on a numerical integration of the two 
equations. These calculations were performed 
on the Model 1103 of Remington-Rand Univac. 

The region of the phase plane which is of 
interest is for 0 < € <1 (i.e. concentration, 
between 0 and 100 per cent of inlet value) and 
1-75 < » < 2-45 (i.e. temperatures between 300 
and 580°K). This region is presented in Figs. 
11-20 and to clarify the way in which the con- 
centration and temperature depend on time a 
number of the trajectories are so plotted in Fig. 10. 

The phase plane when k = 0, shown in Fig. 11, 
has already been the subject of comment in par. 8 


and it is in any case entirely similar to that found 
earlier by AmuNpson and Bitous [2]. The 
trajectory DA going from the bottom left hand 
corner to the low temperature critical point A 
is shown as a function of time in Fig. 10 (a) and 
(b) as the curves labelled A there. The curves 
in the top right hand corner of Fig. 11 which 
cross the high temperature border (such as that 
at E) will go to a high value of » and low value of 
€ before making a sharp turn and coming back 
beneath the separatrix BC to the high tempera- 
ture steady state C. From a practical point of 
view an extreme departure of temperature could 
probably not be tolerated and shows how inadvis- 
able it would be to start the reactor full of the 
reactant mixture at inlet composition. 


Fig. 12 shows the phase plane for k = 0-847 
very shortly before the critical point C becomes 
unstable. It shows little change from Fig. 11 
except in the neighbourhood of the critical 
point C. The approach to C is very, very slow 
the last point of the trajectory as plotted, namely 
D, was reached in time + = 4. The trajectory 
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Fic. 11. Phase plane, k = 0. 


was, however, computed long beyond this to 
7 = 249 and this point is shown £, still not too 
near C. The arm FB of the separatrix has also 
moved over to higher temperatures. At k = 1 
(Fig. 13), the critical point C is definitely unstable 
and the trajectory FGHB shown there winds 
outwards from C to B. The other arm of the 
separatrix BDEA sweeps around very close to 
the arm BHGF but goes up to the stable critical 


0 














iS BO D1 22 23 D4 
Temperature, 7) 


Phase plane k = 0-847. 
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Fic. 13. Phase plane, k = 1. 


point A. Thus the calculations definitely establish 
that the alternative (i) of Fig. 8, which was 
thought probable on topological grounds, does 
in fact obtain. A trajectory from any point such 
as J within the loop will run down to a corner near 
D and then proceed to A passing through the 
very narrow strait between the curves HB and 
DE. 

The change from k = 1 to k = 2, i.e. between 
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Fic. 14. Phase plane, k = 2. 
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Figs. 13 and 14, is comparatively slight and 
entirely as has been predicted. The unstable 
critical point C has moved toward B and its 
instability has become much more obvious. 
The arm of the separatrix BDA still has the same 
form and the trajectory EA is typical of those 
starting within the loop. Attention may be 
called here to the fact that not all trajectories 
have been drawn in full. This is not necessary 
and saves cluttering up the figures, particularly 
later ones. It is clear, for example, that the 
trajectory FG ultimately goes to A on a smooth 
and obvious course. 

In Fig. 15, k = 4°5, at which value the moving 


point C coincides with B, and the limb BC of the 


separatrix vanishes. All trajectories go to the 
stable point A and one such from inside the loop 
is shown, DA. The curves labelled I in Fig. 10 (a) 
and (b) are time plots for the trajectory EA of 
Fig. 15. With k > 4-5 but . 
of the phase plane is similar to Fig. 14 save that 
B (€ = }, » = 2) is now the unstable node and 
C, between B and A, is the saddle-point. At 
approximately k = 5-9 A 
leave behind a stable limit cycle. Fig. 16 shows 
6 just after this limit 


5-9 the appearance 


and C coalesce and 
the phase plane for & : 
cycle has evolved ; note that its top end is still 
very sharp, according with the fact that it is 
formed from a curve with a cusp there. Trajec- 


Concentration, é 


Fic. 15. Phase plane, 


tories from all points outside wind themselves 
around it, while such a one as DEF starting from 
a point within winds outwards and very quickly 
takes up the form of the limit cycle. 

The way in which a limit cycle gets into its 
stride may be seen from the curves of Fig. 10 (c). 
These all apply to the case k = 7, for which a 
phase plane diagram is not given but whose 
limit cycle can be seen in Fig. 21. The letters 
by the curves denote the initial point of the 
trajectory according to the following scheme : 
A,€& € 1-0, » 
C, é ‘ »&€ = 02,79 
E, & = 08, 7 ", € = 055, n 


1-85 


2°25 


2-01 


The fully developed limit cycle is the are WXYW 
to the right of the figure which the trajectory W 
starting from within picks up very quickly. 
Although, as Fig. 21 shows, the part of the limit 
cycle with greatest € is rather pointed when & 
is plotted against + this part is a flat hump. 
This is because the limit cycle is traversed very 
slowly at this low temperature, high concentration 
part whereas at the high temperature end the 
movement is very rapid and produces the sharp, 
almost cusped, dip near the W's of Fig. 10 (ce). 
It is clear that all the other curves are going to 
tend towards the same form, though since they 
start from widely different parts of the phase 
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Fic, 16, Phase plane, k = 6, 
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Concentration, € 


Temperature, 1) 


. Phase plane, k = 8, 


plane there is a phase lag between them, which 
shows up as a displacement of the ultimate form 
of curve in the + direction, 

Fig. 17, the phase plane for k = 8 shows the 
development of the limit cycle. The curve O 
of Figs. 10(a) and (b) shows & and 7 plotted 
separately against + for the particular trajectory 
CGHIJ. A closer study of the development of the 
limit cycle is given in Fig. 21, from which its 
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Fic, 18, Phase plane, k = 9. 
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Fic. 19. Phase plane, k = 10. 


decrease in size and tendency to become more 
elliptical as k +9 is evident. The inset shows 
the way in which the periodic time, 7’, of the limit 
cycle and also its least semi-diameter, B, vary as 
functions of k. It is seen that though the size 
of the limit cycles falls rapidly to zero as k +9 
the periodic time T' tends to a limit of approxi- 
mately 4. This is because the limit cycle is 
closing down on a critical point at which the 
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Fic. 20. Phase plane, k = 30. 
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Concentration, ¢ 


Ternperature, 7 
Fig. 21. Limit eveles. 
velocity vanishes; though the perimeter is decreas- 
ing, the periodic time, which is approximately 
the quotient of the perimeter and velocity, tends 
to a non-zero limit. 
The situation, & - 
Within the stable limit cycle the tendency of a 
trajectory to spiral towards it is almost impercep- 


= 9, is shown in Fig, 18. 


tible and this region is virtually one where small 
oscillations will persist. Even from outside the 
approach is extremely slow as is shown by the 
trajectory CDEFGHIJKL, which even after five 
oscillations is far from the limit cycle. This is 
wholly in keeping with the picture formed from 
the calculations of the preceding paragraph, where 
it was seen that for as low a value of k as 9-025 
the limit cycle system would have vanished. 
Extensive calculations for values of & just greater 
than 9 confirm this. 

The way in which control of the whole phase 
plane is finally attained is seen in Figs. 19 and 20, 
In the former, for k = 10 a rather weak form of 
control obtains, and its takes a number of oscilla- 
tions before the concentration and temperature 
settle down to their controlled values; the curves 
Q of Fig. 10 (a) and (b) show this. In the latter 
when k = 30 all trajectories spiral in to the 
controlled point and their oscillations are almost 
immediately negligible. This is confirmed by the 
curves T of Fig. 10 (a) and (b). 


CONCLUSIONS 


It has been shown how the heat transfer to a 
cooling system subject to control can be approxi- 
mated by a quadratic expression in the reactor 
temperature. Using this approximation a simple 
system with a dominant first order reaction is set 
up for study by the topological methods of non- 
linear mechanics. The uncontrolled system has 
three steady states, two stable and a third 
unstable, of which the control is being studied. It 
is shown that as the control parameter is increased 
the steady states change both in position and 
stability, and that stable limit cycles can develop. 
This represents a particularly undesirable feature 
of reactor behaviour corresponding to a con- 
tinuous cyclic variation of temperature and 
composition, which the reactor ultimately follows 
whatever its starting conditions. A more subtle 
difficulty arises just after local control of the steady 
state is attained, for this stable limit cycle, 
Instead an 
unstable limit cycle is generated around the 
steady state and within the stable limit cycle. 
Thus, although the criteria of local stability assure 


though small, does not disappear. 


control at the steady state, this control only 
extends to small disturbances which do not take 
the system outside the unstable limit cycle. 
This phenomenon emphasizes the importance of 
using the method of non-linear mechanics, 
without which it would go undetected. Moreover, 
the fact that it happens even with a perfect 
controller (though somewhat ephemerally in the 
particular example considered) suggests that it is 
still more likely with imperfect control. 


List or SYMBOLS 
Heat of reaction (2) 
Temperature 
VU* /qaq(— Al) reduced h.t.f. (2) 
Volume of reactor (2) 
Chemical species (2 and 14) 
Abbreviation for the function 
1+ 1 B 


io. 
exp fe (; _ :)| (9) 
(1 + «)(1 + B) — y (3) 


~ ¥ (3) 
Proportionality constant 
: Flow rate of reactants (2) 
Ye = Flow rate of coolant (2) 
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t = Time (2) r = qt/V, reduced time (2) 
P = Reduced reaction rate (2) w = Roots of characteristic equation (3 and 6) 
a, 8B, y = Stability parameters (3) Suffixes and affixes 
7 = epT /1%g(— All), reduced temperature (2) c¢ = Coolant 
"y = kK"! %es/q, half range of temperature causing 0 = Initial value 
valve to move 8 = Steady state value 
é=2 Io reduced concentration (2) * = Deviation from steady state (4) 
A = (k — k*) in paragraph 10 m = Maximum value 
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The principles of programming reactor calculations. Some extensions 
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Abstract—The increasing use of automatic computers for chemical engineering calculations 
makes it desirable to give an outline of the programme used here. This is so constructed that it can 
To round off the question of proportional 
Finally the 


extension is given to simultaneous reactions in more than one reactor and the general conclusions 


be rapidly modified for many more difficult problems. 
control with very high amplification the limiting form of off—on control is discussed. 


of the problem are discussed. 


Résumé L utilisation croissante des calculateurs automatiques dans les caleuls de Génie 
chimique rend nécessaire de donner un apercu du programme employe ici. Tl est établi de fagon a 
pouvoir étre modifié pour les problémes difficiles. Pour circonvenir la question du controéle 
proportionnel avec une trés grande amplification, les limites du contréle “ tout ou rien” sont 
précisées. Finalement les auteurs donnent lextension de la méthode au cas des réactions simultanées 


prenant place dans plusieurs réeacteurs et discutent les conclusions générales du probléme. 


Zusammenfassung Die steivende Verwendung automatischer Rechner fiir verfahrenstech- 


nische Berechnungen lisst es wiinschenswert erscheinen, das hier verwendete Programm zu 
beschreiben. Dieses ist so gestaltet, dass es schnell fiir die verschiedensten Probleme abgeindert 
werden kann. Um die Frage der Proportionalregelung mit sehr hoher Verstirkunyg abzurunden, 
ist die Grenzform der Aus-Ein-Regelung behandelt. Schiiesslich werden die Betrachtungen auf 


gleichzeitige Reaktionen in mehr als cinem Reaktor ausgedelhnt und die allgemeinen Folgerungen 


des Problems diskutiert. 


12. Tue PrRIncrireLes OF PROGRAMMING 


rue CaLCuLATIONS 


Tue considerable amount of calculation which 
lies behind the results presented here was per- 
formed on the Univac Scientific Computer, 
Model 1103. This machine is well-suited, both 
in size and versatility, to this type of calculation, 
which would indeed have demanded an impossible 
amount of time and labour to do by hand, 
While much of the detail of the programme is 
particular to the machine used, and thus un- 
suitable to be described here, its overall strategy 
and general principles would apply to almost 
all modern high-speed computing devices and 
are the same for the more general problem 
treated later. In view of the increasing importance 


of automatic computers in this field we sketch 
the programme in such a form as to allow it to 
be adapted to any problem of performance or 
control of a continuous flow reactor. 


the concentration 
of the 
integrate the pair of equations (9.1) and (9.2) 
from a given initial state €,, », until such time 
This 


may be either because the system is very close 


To obtain temperature 


histories converter it is necessary to 


as the future behaviour is predictable. 


to a steady state or because it is performing a 
limit cycle. The Runge-Kutta method of integrat- 
ing non-linear simultaneous equations numerically 
is well-known (e.g. Ince [7]+) and has been adapted 
to machine use by Giix [6]. A finite interval 
Ar is chosen and the corresponding increments 


*For ease of reference the numbering of the sections and figures is continuous with Parts I and II. 


**Present address : 


The Mathematical Institute, University of Edinburgh. 


TAIl references are listed at the end of lVart I of this paper (p. 131). 
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A€é and Ap» are calculated from the equations in 
such a way that the error is proportional to 
(Ar)®; since the number of steps in a trajectory 
is inversely proportional to Ar the error in the 
whole trajectory is of the order of (Ar)*. 

As each point of the trajectory is obtained 
several tests are applied to decide whether (a) 
a monitoring signal should be printed out on the 
typewriter, (b) the calculated point should be 
permanently recorded or (c) the trajectory should 
be stopped. 

The first class of test is merely a convenience 
to give the operator some idea of how the trajec- 
tory is shaping; it was arranged for a period to 
be printed for every 25 points of the trajectory 
and for the sign of » to be printed every time it 
changed. In general far more points are calculated 
than can possibly be used and since the recording 
of a number involves converting it to decimal 
and punching a tape, both relatively lengthy 
operations, it is uneconomical to record more 
than is necessary. When a point £», yp had been 
punched out it was therefore arranged that no 
point would be recorded until lé Ep| or 
ly — np| exceeded certain limits. Finally the 
stopping tests can be worked in with the preceding 
tests quite easily. If the number of points between 
each punch out is recorded then the growth of 
this number beyond a certain limit M denotes 
that the steps are so small that the movement is 
almost negligible and the trajectory is approach- 
ing a critical point. Again the change of sign of 
7 from + to denotes the fact that » has 
passed through its maximum value, ie. the 
trajectory has turned a corner. If this maximum 
value is compared with the value of » the previous 
time that it passed through its maximum the 
difference will become very small as a limit cycle 
or focus is approached. The trajectory may be 
stopped therefore when this difference falls below 
a preassigned value. A trajectory is also stopped 
when it goes outside a prescribed region of the 
phase plane. In all cases the last point is recorded 
and a monitoring sign to show the type of stop is 
printed out. The details of these tests will be 
published elsewhere. 

The programme is composed of a number of 
sub-routines and these are joined together in the 


scheme shown in Fig. 22. There are three routines 
which from given values of € and » will perform 
the necessary substitution in the equations : 


A calculates P (£, ») from a given kinetic 
expression, 

B calculates € from equation (9.1), 

C calculates 7 from equation (9.2). 


[Operator entry | 


a |Jo 


prey) 


General driving 
routine 


+4 
307 30<-70 


*—"_ Calculates 
Initial | 


fandn 
|Applies tests 


Calls for 





Fic. 22. The programme. 

These three are organized by a driver sub- 
routine D, which calls for them in the correct 
order. The core of the programme is the general 
driving routine E, which calls for the calculation 
of € and » and then for the next point of the 
trajectory to be calculated by the sub-routine 
G; G is the sub-routine for carrying out GrILL’s 
adaptation of the Runge-Kutta method, in the 
course of which it calls on D for other values of 
the derivatives. The general driving routine E 
also applies the various tests described above and 
accordingly instructs the following sub-routines : 

H, the stopping routine arranges the necessary 


print out and then calls for the next 
trajectory, 

I, the typewriter routine prints out the monitor- 
ing symbols and the indication of the type 
of stop, 

J, the punch-out routine, which converts the 


given point to decimal and punches a tape. 
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Finally there is a routine F which is called for 
at the beginning of each trajectory either by the 
operator or by the stopping routine H. This deals 
with the setting up of certain registers to standard 
values and inserting the initial point of the 
trajectory. Fig. 22 shows the interrelation of 
these ; F EH is the main loop of the control which 
is traversed for every trajectory. 

This scheme has been given in some detail for 
it will apply to most types of computer and its 
reduction to these basic units of organization 
allows the following generalizations to be made. 
If the kinetic equation of the reaction is changed 
it is only necessary to alter the sub-routine A. 
If in the equations one of the quantities which 
has been taken as constant (e.g. g) is made time 
dependent this only calls for a further sub-routine 
to calculate it and a slight modilication of D to 
organize it in relation to B and C. If more than 
one reaction or reactor is considered further 
sub-routines for calculating the derivatives must 
be constructed and D modified to organize them. 
The sub-routine G can be written for any number 
of simultaneous first order equations so that no 
extra difliculty is introduced here, the tests used 
in E will, however, require to be made more 
general, 

Finally from an economic standpoint the use of 
a large scale computer is easily justified. One 
trajectory involved 1320 steps and took 200 sec 
including the punching out of some 260 points. 
A skilled human computer working fast might 
perhaps get 6-6 steps done in Lhr, ic. take 
3600 times longer. Machine time rented at 
$360/hr is thus equivalent to human beings at 
10 ¢/hr. 


13. Orr-on CONTROL 


In par. 7 it was seen that as the gain k was 
increased, the range of temperature over which 
the flow rate was controlled decreased to zero, 
so that in the limit we had an off-on control. 
The equations in this case may be written 


é=1— £¢— P(é,n) 
7 = — 7 — Plé.4) + hy, — 9) 


where 


(18.1) 
(18.2) 


0 1<% 
k, 7 = %, 
kg 17>, 
If £,, n, is to be a steady state, k, must be such 
that 


Yo ~ % P(E,, s) 7 k(n, cam "s) =0 (13.4) 


Further, if the phase plane diagram for k = k, 
is drawn it is both necessary and sufficient that 
the half plane » > », should contain no stable 
critical point for then trajectory from 
this part of the plane will cross the line 7 = », 
towards But the 
behaviour of the system for » < », is represented 
by the phase plane for k - 0 it is again necessary 
and suflicient that this should contain no stable 


hn (18.8) 


any 


lower temperatures. since 


critical points for » < »,. 
Solving é = 0 for £ = £(») and substituting 
in» =O 
P(E{n), 9) 
= J(n), say, 


(1 +) 7 Ve —kn, 
(13.5) 
of which the right hand side is a sigmoid curve as 
shown in Fig. 23. The left hand side gives three 
straight lines for the three values of k, namely : 
for k 
for k 

through », = 
for k = k, 


through », 


0 a line (A) of slope 1 passing through »,; 
=k, a line (B) of slope 1+, passing 
(no + By n¢)/(1 + &); 
a line (C) of slope 1+, passing 
(no + Kz n-)/(1 + ky) 


Consider also the lines: D, that tangent to 
the curve F of slope 1 whose point of tangency 
is of least y, and E, the tangent to F of slope 


4 








{ —T 2 
c M5", N 


Fic, 23. Off-on control. 
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1+, with greatest value of » at tangency. 
Let these intersect the axis f = 0 at n, and 7, 
respectively. Then provided , > 4, the line A 
must intersect F once only and that at a tempera- 
ture greater than »,. The line B intersects F at 
the point » = »,, and again providing », < n, 
the line C will intersect F only at a temperature 
lower than »,. These three conditions provide 
suflicient restrictions on the disposable parameters 
to ensure only one steady state in the phase 
plane. They are 


(n, - No) + k, (n, = Ne) = S(n,) 
No > 14 
to + he ne < (lL + he) m2 


Since k, will generally be fixed by the limit of 
flow the last condition gives an upper bound for 
7. the coolant temperature and the second a lower 
bound for m,, the inlet reactant temperature. 
It is clear from the geometry of Fig. 23 that while 
these are sufficient they are not always necessary. 
Off-on control has been treated as the limiting 
case ; in practice it would be a chattering type 
of control which would be unpleasant but the 
conditions (13.6) may be applied to the function 
U to insure stability even with the highest gains. 


(18-6) 


Tue E@vuations ror Many ReEaActTIoNs 
AND REACTORS 


14. 


Finally the most general situation with a chain of 
stirred tank reactors will be discussed with the intent to 
show what equations need to be set up to describe the 
behaviour of the system. As has been remarked the general 
principles of this analysis apply to the general situation 
and some features, such as the programming of the calcula- 
tion, require only trivial modifications, The phase plane 
presentation is, of course, lost as soon as more than two 
variables are considered and the corresponding representa- 
tions in space of higher dimensions requires a little more 
imagination. 

Consider first the situation when the inlet concentrations, 
temperature and flow rate are time dependent. A mass 
balance for each substance gives the equations derived 
in par. 2, 


Wty 2) - VR (14.1) 


(14.2) 


It would still be possible to eliminate all the concentrations 
except one, namely x. For multiplying (14.1) by A, and 
adding the two equations gives an immediately integrable 
first order equation 


t 
—fava TF 
° [ @ 
V 


0 


A faoft) r(t)) 


T 


1iryit) — Toft) + riglt) — Fol a} {« xp |q ‘Vdz)dr (14.3) 


0 
the necessary modification of the relation (2.2) to which it 
In this case the 
taking 


reduces if apt) and xjo/t) are constants. 


variation in flow rate is 


t 


accommodated — by 


T q(t) dt /V though the reaction rate P then becomes 


time dependent. The elimination (14.3) is 


unprofitable from a practical point of view, for this formula 


given by 


would in any case be computed by a process equivalent 
to the integration of the differential equation. However 
because the elimination is possible there are essentially 
only three variables, one concentration, the temperature 
and the time. After a reduction to dimensionless concentra- 
tion and temperature analogous to that of par. 2, there 


are two equations of the form 


dé ~ dy » 
7 P(£, 9. 7), - QE, 9, 7) 
dr dr 


These could be written in a more symmetrical form by 
the device of replacing + on the right hand side of the 
equations by a dependent variable {, and adding a further 

dt 


rT or — 1. 
dr 


equation { = Then, using the dot notation 
for d/dr 
Wéi.n0. © = RE 


The behaviour of the reactor is 


é P(£,9,0). 9 


where of course R = 1. 
thus represented in a three-dimensional space, through any 
point of which one trajectory passes. A projection of these 
trajectories on the £, plane will retain some similarity 
to the time-independent phase plane, but these projections 
will intersect at many points besides the critical points. 
If the time dependent quantities tend ultimately to steady 
values then the trajectories will all be asymptotic to lines 
parallel to the { axis through the critical points of the 
phase plane which correspond to these steady values. 
In particular, if trajectories are taken backwards from the 
points of the ultimate separatrix they will intersect the 
plane { = + = 0 ina curve, which will separate the regions 
that are attracted to one or other of the steady states. 

It will now be shown that when m independent reactions 
are taking place simultaneously in a stirred-tank reactor, 
the m mass balances and the heat balance are sufficient to 
describe the system. Let the number of chemical species 
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taking part be n (>> m) and denote these by X;,j = 1,2,.. 
n. If the stoicheiometric coefficient of X; in the ith reaction 
is Aj;, the equation for that reaction may be written 
n 
2 Ay; Xj = 0. 
j=l 


independent there is no linear combinacion which expresses 


Since the reactions are assumed to be 


any one of them in terms of the others and the species 
may be reordered so that the first m equations must 
include at least one of X,,. . X,,. 
of species the equations may then be reordered so that 
no A 


With this reordering 


a= 0, i=1,...m. The equations retain their 


meaning when divided by a constant so let the ith equation 
be divided by Aj; and write aj; = Ajy/ Ay (ie. ay = 1). 
The equations may now be written in a standard form 


O,7 1,...m, or using the matrix notation 


* Qn xX, 


++ 4m a m+ ** 


42n X, 


coalietliaiis eee 
=[4\ B)x=0 











oat Om = 1 6 Gymiis + 4mm | Xn (14.4) 


@ is the m x n matrix of stoicheiometric coeflicients aj; and 
A isthe leading square m x m submatrix, which by arrange- 
ment is non-singular and has units along the diagonal. 

If N; moles/unit volume of X; are formed in the ith 
reaction then by the stoicheiometry of this reaction aj; N; 


moles of x; are formed and hence the total number of 


m 

moles of X; formed in all reactions is 2 aj N,. Let 2; and 
t= 1 

instantaneous initial concentrations 


Ty be the and 


(moles/unit volume) of X; and write 


— 


7) 


and N 








Then clearly 2; = rj yj Nj OF in matrix notation 
tel 
x = Q’N, where a prime denotes the transpose of @. 


Using the partitioned form of @ and z 

ye An | 
(14.5) 
z=BN | 


Since A is non-singular we can find N in terms of y and 
subsitute in the equation for z obtaining 


z=B Ay (14.6) 


This expresses the concentrations of X,,,,... 
linear combination of those of X,,... - Xin: 

If (dx;/dt); denotes the rate of change of concentration 
of X; due to the ith reaction alone, the rate of this reaction 
can be uniquely defined by a rate equation of the form 


X,, as a 


1 di; 
~ | “)) = RAT, 24... .%_) = RAT. 24, -++%—_) (14.7) 
aj dt i 
i eS 


A matter balance for xX; gives 


- : 
— Tyo) + = ay R;. j 


It is now claimed that only first m of these equations are 
independent for one may write 


(14.9) 


(14.10) 


Then substituting from (14.6) into (14.10) 


(5, ‘y BA“ y — BR 


at 


BA (i + r)y . an =0 


and this equation is identically satisfied in virtue of 
(14.9). 

If q and jg are functions of time one may multiply 
(14.9) by B’A" and subtract (14.10). Then 


d q d q — 
(i+ pr)? (a + p) (aA ¥) 


may be integrated to give z in a form similar to (14.3). 
Again, however, this is only necessary from a theoretical 
point of view to establish the dependence of the last 
(mn — m) equations on the first m, for computation all n 
equations would be used. 

If h; is the heat released in the ith reaction by the 
formation of one mole of X; then a heat balance gives 


m 
-T)+V SAR, -— UV 


i-1 


_ a 
Vee—= = qep (Ty (14.11) 
where U* is the rate of heat removal by the cooling 


system. A reduction to dimensionless concentrations, 
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£;, and temperature, », may be effected entirely similar 


to that of par. 2 and obtain 


) hd > 
+ = ayP; 
| 


sj 
(14.12) 
= b;P; 
i 
where the constants 6; are proportional to the heats of 


reaction /h;. If ¢ {£y’- ++ Em Mer Co , are 
1 


" ~ U} and Cc 5 


primes denoting transposition, then 


| P P(t) = {P,,... P,,}, the 


(14.13) 


Though these equations are written in matrix notation 
it must not be thought that they are linear. In general 
P will be anything but linear in the variables (. The 
steady state values are again found by solving the non- 
CP = 0. 
(, then 


linear equations ( ra If ¢, is a solution of 


- 


these and 6 = ¢ 


0 (8*) 


grad P - 


(14.14) 


al 4 
» 


and 0 (8) represents terms of the order of the squares of 
the differences { — [,. The stability of the steady state 
is thus entirely defined by the latent roots of the matrix 
D — I. If one writes 


aP; 
a = ry; , general parameter corresponding to a of 
“*) /apar. 3 


\ 


aU 
- |—] , the £ of par. 3 
i”), 


d 


-- ) , the parameter corresponding to y of par. 3 
s 


| 


The stability of the steady state will thus depend on at 
most (m2? + m 1) parameters. 


find the conditions under which perfect 


It is necessary now to 
proportional 
control on deviations of temperature will always succeed 
in locally stabilizing a steady state. 

Perfect proportional control of the heat transfer function 
on signals of the deviation of temperature is equivalent to 
modifying the value of 8. For if q, is no longer a constant 
in the function U (», qg, n-) then 


(-) (~ dq, 
On, & i dy 
If the latent roots 
local stability is 


and dq,/dn = k, the control constant. 
of (D 1) 


assured. Let 


have negative real part 


I D M N where M 


“Cc [xy - vi] 


- Pp 


and denote the leading m » 
expansion of Al+ M+N 
(m + 1)th degree in A. 


m submatrix of M by L. The 
0 is an equation of the 


Am + +” sp, (M + N) + A*™—! spg(M + N) 


MsN =0 


where sp, A is the sum of all the diagonal submatrices of 
order r X r of a matrix A, 


det (M + N) 
M+N M' 


SPm a (M + N) 
BL 


Similarly 


spy (M + N) spyM B, sp,(M + N) 


= sp,M + 8B sp, “f L 


Writing the characteristic equation as 


ym 


+ a, A® + ag rA* Pe sil 0 


m+1 


The Routh-Hurwitz criteria for roots of negative real 
part are 


ay > 


and from the rth criterion 


1 O O 
ag 1 
as Ug deg 


there emerges 8’ multiplied by a factor (spyL) (spgL).... 


(sp,_,L). This can be made to dominate the determinant 
and force it to become positive provided all the sp,L... 
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This sufficient condition involves 


 —— I" + a 
Gan, > +: 1 —_ . 


Sp, are positive. 
only the aj for 


L = (yj) 
im 


where |; 2D ay; ap. It is interesting to note that the 
kel 


possibility of successful control depends only on the me 
parameters 2,;, but the ease of control will depend largely 
on the further m parameters y;. 

If r reactors are connected in series so that the product 
of one passes (perhaps after modification) to the next 
there will be a set of (m 1) equations for each reactor. 
In addition relations between the inlet concentrations and 
temperature of one reactor in terms of the exit conditions 
of the preceding one will exist. If an aflix (4) denotes the 
value appropriate to the Ath reactor, a straight through 
connection with no heat loss between one reactor and the 
next simply gives 2{**) = ay, T&-) — 7. It 
will be supposed however that in general some of the 
flow is bypassed and the flow rate made up from some other 
stream. Let q® be the fraction of the flow q*) out of 
the (k)th reactor which goes to the (4 1)th and let an 
of tlow (*) and 


a be introduced at this point. Then if mixing is complete 


additional stream rate q composition 


(k+1) 
0 


a: (k) @ A) gi*t+)) 


q 
and 
q* +1) 


q*®) a (*) 


q 
If the concentrations 2 are reduced to dimensionless form 
é then 


tig"* +1) = i*) &*) «il xi*)) é{*) (1 4.15) 


and the «) = 1 for a direct connection. In a similar 
way the reduced temperatures will be connected by a 
relation of the form 


no*t® 7 \) f*) al fd) 7 (*) ul*®) (14.16) 


where «) is a heat loss between the reactors. 
that 
necessary and sufficient to describe the behaviour of the 


These equations show r(m 1) equations are 
chain of r reactors in which m reactions are simultaneously 


taking place. If however the inlet concentrations and 
temperatures are time dependent it is not wise to perform 
the possible eliminations but use the full set of 
(2rn + r — n) equations of the form (14.12), (14.15) and 


(14.6). 
CONCLUSIONS 
An outline has been given of the strategy of 
the programme needed by a computer of the 
modern high speed type for calculating the 
transient behaviour of a continuous flow reactor. 


An important feature of economical programming 
is a flexible system of stopping tests, and one such 
is described. Simplifying assumptions have only 
been introduced for the purposes of exposition 
and the final paragraphs give the way in which 
certain extensions may be made. With very 
high amplifications of the controlling signal the 
controller meets the limits of its movement almost 
instantaneously and to ensure stability in the 
whole phase plane the necessary conditions for 
stable off-on control are discussed. Finally the 
method of extension of simultaneous reactions 
and many reactors is outlined, and it is shown 
that similar stability criteria can be established. 
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Abstract—A method is proposed whereby the flow structure of a system, or the physical structure 
of a packed bed or porous solid, may be analysed in terms of a model chosen (a) to approximate 
hydrodynamically to the system under investigation and (b) so that its effeet on a flow of fluid 
containing a sinusoidally varying concentration of solute can be mathematically analysed. 

The attenuation and phase angle of the output stream are measured for a number of different 
frequencies ; as many values of the frequency being chosen as there are values of the unknown 
variable to be found. From the measured results linear simultaneous equaticns are set up, which 
can be solved to give a histogram of the distribution of some desired parameter. 

Results are given for models which illustrate the method and which are likely to be of use in 

practice. 
Résumé—L‘auteur propose une méthode qui permet d’analyser le mode d’écoulement d'un 
systéme, ou la structure physique d'un garnissage ou d'un solide poreux, en fonction d'un modéle 
choisi (a) pour avoir une bonne approximation hydrodynamique du systéme étudié et (b) pour 
analyser mathématiquement les effets produits sur écoulement d'un fluide contenant un soluté 
dont la concentration varie de facon sinusoidale. 

La diminution de l'amplitude et le déphasage du courant sortant sont mesurés pour 
diverses fréquences. On choisit autant de valeurs de fréquences qu'il vy a de valeurs inconnues 
a trouver. Partant des résultats linéaires experimentaux, on établit des équations simultanées : 
équations qui peuvant étre résolues pour donner la variation en fonction du temps de n’importe 
quel parameétre désireé. 

Les résultats sont donnés pour des modéles qui illustrent la méthode et qui sont vraisem- 
blablement utilisables dans la pratique. 


Zusammenfassung—Um das Stromungsverhalten eines Systems oder die physikalische 
Struktur cines Festbettes oder eines pordsen Korpers zu erforschen, wird die Verwendung eines 
Modells vorgeschlagen, mit dem man (a) das betrachtete System hydrodynamisch annihern kann 
und das es (b) méglich macht, sein Verhalten bei sinusfOrmiger Veriinderung der Eintritts- 
konzentration analytisch vorauszuberechnen. 

Die Verdinnung und der Phasenwinkel des austretenden Stromes werden fiir verschiedene 
Frequenzen gemessen und zwar werden ebenso viele Frequenzwerte gewiihlt wie Werte der 
unbekannten Variablen zu finden sind. Aus den Messergebnissen werden lineare simultane 
Gleichungen aufgestellt, deren Losungen ein Spektrum der Verteilung eines gesuchten Parameters 
ergeben. 

Zur Darstellung der Methode werden die Ergebnisse fiir Modelle mitgeteilt, die fir den 
praktischen Gebrauch geeignet sind. 


INTRODUCTION of this flow-structure, i.e. of the distribution of 
Ix many systems, in the field of chemical engineer- residence times, is often needed. For example, 
ing and in others, there is flow of fluid in which the flow-structure in a reactor is needed for 
the bulk flow is made up of the sum of a number the analysis of its behaviour and for its design. 
of individual flows, so that the “ residence-time ” Again, considering the system itself, e.g. a 
of a small portion of the fluid may have any porous solid or packed bed, this is often treated 
value between zero and infinity. A knowledge as though it consisted of a collection of pores of 
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various lengths and diameters. To get a complete 
picture of this model it would be necessary to 
determine the frequency distributions of the 
lengths and of the radii of the pores. Instead, 
one has had to use an average value of the lengths 
of the pores (obtained by multiplying the length 
of the bed by a tortuosity factor), and an average 
value of the hydraulic radius m, obtained from 


€ S, 


n= 


where ¢€ = porosity, 


S = specific surface. 


Even this average value is not of use when there 
is a wide spread in the distribution of pore 
diameters [1]. 

In this paper methods are suggested whereby 
the frequency distributions of lengths and of 
of the 
considered to be the hydrodynamic equivalent of 
the packed bed or porous solid may be found 
experimentally. 


diameters pores in a model which is 


In addition, the problem of residence times is 
considered for the special case where the fluid 
flow is considered to occur in a number of more 
or less identical channels ; communicating with 
each channel are a number of stagnant pockets 
of fluid. The liquid in these pockets has infinite 
residence time but would play a part in any 
chemical or physical process by the diffusion of 
material into or out of them. 

A method is suggested whereby the distribution 
(by volume) of pockets of different depths, and 
hence the amount of “dead” liquid relative to 
the total empty space in the system, may be 
found, 

These are two examples of the type of informa- 
tion which might be required, and in the body 
of the paper will be developed briefly the 
mathematical equations which apply to these 
specific models. The same type of approach 
can be used for other models, the choice of which 
would have to be made after consideration of the 
actual system to be studied, and the procedure is 
summarized in the section headed Conclusion. 

The problem is made up of two parts; the 
first is to determine times and the 
second is to relate this knowledge of residence 


residence 


times to the flow structure within, or the physical 
structure of, the system. 

Danckwerts [2] has considered the use of 
step-functions and impulsive or “delta-functions” 
(Dirac) for the purpose of determining residence 
times but in this paper the use of “ frequency 
response will be considered because there are 
additional unknown variables in these problems, 
(each point on the distribution curve is an 
* unknown ” which it is desired to find), and each 
additional value of the frequency used enables 
another equation to be set up. 

In recent years a number of papers have 
appeared dealing with the application of frequency 
response to the analysis of flow systems. Thus 
Rosex and Winscue [3] have dealt with the 
theory of its application to chromatographic 
columns; Detster and Witwerm [4] with a 
theoretical and practical study of its application 
to a bed of porous solids ; McHenry and WILHELM 
[5] with its use in measuring mixing in a bed of 
spheres; Kramers and ALBerpa [6] with its 
application in continuous flow systems and Keyes 
[7] uses it in measuring radial mixing in turbulent 
flow in a tube. 

To deal with the second part of the problem, 
viz. the relating of the residence times to the 
flow structure or physical structure of the system, 
the effect of any chosen model on a sinusoidally 
varying property must be calculated. Here this 
property will be the concentration of some solute, 


and the results for some models will be presented. 


Tue Errecr or THE 
PHASE 
OF THE 
Model 1. 
The structure is a channel or pore having a 
uniform cross-sectional shape and area; com- 
municating with it are dead spaces or “pockets” 
distributed uniformly along the length of these 
channels. There is no flow in or out of these 
pockets and any transfer of solute in and out 
is by molecular diffusion only. 

Here, then, the total volume occupied by 
fluid is made up of a part, in the channel, of 
which the residence time is finite, and part 
(the pockets), in which the residence time is 


MODELS ON THE 
AMPLITUDE-MODULATION 
CONCENTRATION WAVE 


AND 


Channel with dead-spaces 
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infinite. It is required to find the relative amounts 
of these two parts. 

The method of attack is based on that of 
Deister and Witne to. Fig. 1 shows a channel, 
distributed uniformly along which are pockets 
all of length / normal to the axis of the channel. 


Pocketswolume -& 
c/S crea: jy’ 
1. Model 1: a channel with distributed pockets 


of stagnant fluid 


The assumptions made are that : 


The channel and all the pockets are full of the 
one fluid, 

The pockets are of prismatic cross-section, 

The concentration of solute is uniform across 
the channel, 

The width of the pockets, in the direction in 
which z is measured, is small relative to the 
wavelength. 


The “pockets”, which have to be given 
geometric shape for the purpose of the analysis, 
mean in reality “ dead-water” spaces in which 
there is no movement of liquid. The results of 
any calculations will therefore be in terms of 
the equivalent pockets. 

At a distance z from the entrance to a channel, 
measured along its length, at time t, the equation 
of continuity for the dissolved solute is 
p 2X, _UdX, _ 2X, 


c 
2? dz at 


—p4. 0» (1) 


where D = longitudinal diffusivity, q = average 
concentration in the pocket (mass volume), 
X, = excess of concentration of solute at distance 
z over the time - average value for the system, 
U = average velocity in the channel, « = volume 
of pockets per unit length of channel, a = cross- 
sectional area of the channel and ff = «/a. 

The solution to (1) is to be given by the real 
part of 


X, A, exp [ i (wt - ¢.)], 


where ¢, = phase angle in the channel at distance 
z from the entrance. The general solution may be 
written 

(2) 


where A, is a complex amplitude. 

In order to solve (1) explicitly it is necessary 
to obtain an expression for dg, dt, and this in 
turn requires that the differential equation 


DY X 


ry. © é 
“y ty (3) 


A dy? 


(X,, is the 
excess concentration of the solute at a distance 
y from the bottom of the pocket over the time- 
average value, and DP is the molecular diffusivity 
of the solute.) 

The solution to (3) will be of the form 


should be solved for the pockets. 


NX. = A, P exp [ 


“¥ 


i(wt — 0,)], (4) 


where the real part only is to be considered. 
@, is the phase-angle in the pocket at a distance 
y from the closed end, and P is an attenuation 
factor relating the amplitude in the channel at 
distance z to that in a pocket, at a distance y 
from the closed end, which is situated distance 
z from the entrance of the channel, as shown in 
Fig. 2. 





AO 














Entronce 


2. The locations of amplitudes and phase-angles 


in a channel and a pocket. 


Equation (3) is subject to the boundary con- 
ditions that at y= Il, X,, = A, cos(wt — ¢,) 
and at y= 0, dX,,/dy = 0 (since there is no 
diffusion through the bottom of the pocket). 

A solution to (3) obeying these boundary 
conditions is obtained by substituting 
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‘ , feosh w'l cos w'l sinh wy sin wy 
an —$ $$ 
leosh w'y cos w’y cosh w'l cos w'l 
cosh wy cos w'y sinh wl sin wl) 


sinh wy sin w’y sin w'l sinh wl) 


, ” 
+ cos 2w'y\'/* 


cosh 2w'y 
and P | 


cosh 2w'l cos 2w'l 


into equation (4), 


} P w 
where w ° 
wre ¢ 2D 


Now °% - D 2S») 


ot i OY Jy! 


vx, :' — 
and | ~ | may be obtained by differentiating 
OY Jy=t 


(4). 


In this way we obtain eventually 


7. A.(Dw'/1)(Q + iR)exp[— (ot —4)], (5) 
j 2w'l — sin 2w’ 
en sinh 2w'l — sin 2w'l 


cos 2w'l + cosh 2w'l 


sinh 2w’l sin 2w'l 


and R 


cos 2w'l + cosh 2w'l 


. ap De Dake 
The quantities ?, — end 

dz oz* 

found by differentiating (2). 
Substituting these into (1), and combining 


with (5), the equation 


PA. 


dz* 


may be 


U dA, 
D dz 
, B ) - 

m4 |. J i(* - *) A.=0 (6) 

D l w' l . 
is obtained for all ¢ after the steady cyclic state is 
reached. A solution to (6), expressing A, as a 
function of z, and involving two arbitrary 
constants, may quickly be obtained. Remember- 
ing that 

X, = A, exp (—iwt) = A, exp (i¢,) exp (—iet), 
and that the boundary conditions 

X, = Ay exp ( a iwt) 
and 

X,->0 


must be satisfied, we get 


x, = Ayexp [up — /(8+ 


 «([(S—F . 
seed mail J ; -)| a (7) 


U? Bw 


im * p| 


where 


F = “)> (8) 


» Bw R 
a= 5 — 80/8), o 


(10) 


o = 2’ l, 


S =4/ (F? + G). (11) 
(One of the arbitrary constants is proved equal 
to zero by making use of the fact that 
U/2D <4/[(S +- F)/2]}.) 

The variations of Q/o and R/o with o are shown 
in Fig. 3. By expanding each in the form of a 
series it can be seen that, to within 4 per cent 
if ¢ < land to within 1 per cent if ¢ < 0-7, Q/o 

-o°/6 and R/o = 1. On the other hand, if 
o> 8, Q/o = R/o=1 











3 4 » © 
oO= 2Qwil 


Fic. 3. The variation of R/a and Qe with e. 


To proceed to the general case ; if we assume 
that, associated with unit length of each channel, 
there is not volume « of pocket of depth, J, but a 
volume «, of pockets each of length 1, as shown 
in Fig. 4., 
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Unit volume of channel 





Fic. 4. A channel with relative volume 8, of pockets of 


length J. 


then, by a slight extension of the above reasoning, 
we have 


xX, - Ay exp | 12D - JF 
alia JOS 


i(wt — ¢.)]. 


p (12) 





A, exp [ 
where 


Uw B, 
4D? D 


D (1 < po oy * 


4/((F*? + (6*)], 


F* = Q,. 


2w'l, 


G* = 


(16) 


and the problem is to determine the value of 8, 
for given values of /.. 

The method of analysis. 
there are a known number of identical channels 


It is assumed that 


in the system of length L, each with its associated 


pockets. Then from (12), putting z — L, 
41) 


a) ls JF 
= /(F}* : “i 


Both quantities on the L.H.S. are measurable, 
so that the values of F*, S*, and hence of G*, 
can be calculated, (assuming that U and D are 
known). 

Hence from equations (13) and (14) the values 
of 28, (Q,/2w'l,) and 2 B,(R, 2w'l,) may be 
calculated. Denoting Q;,20'l, by y 


In (17) 


(18) 


and 


jr 


2 B(Q;, 2w,'l,) by yj, (where j signifies that a 
quantity depends upon the value of the angular 
velocity w;), the steps in the calculation are to 
choose n values of L.. 

From the measured values of the attenuation 
and phase-angle shift, at any value of w,, the 
value of w;, and the values of »,;, for all values 
of r, are calculated. 

This procedure is repeated for other values of 


the frequency, w,...@;...@,. From these a 


j n 


set of n linear equations may be written down ; 


B,, Nin 
B, Nin 


My Bi m1 (19.1) 


(19.j) 


B, Yan 


and the values of 8, corresponding to L, for r 
having the values 1 to n, may be found by setting 
up the determinants of the coefficients in these 
simultaneous equations in the usual way. 


Bn 3, Nnl (19.n) 


* dead- 


2’ B,). 


Further, the fractional volume of 
space " to total empty space is 2 8,/(1 
Model 2. 


diameters 


The model consists of a number of channels 
of different lengths and radii. 
pockets of stagnant liquid. 

Let there be n,, channels of length 1, 
radius R,, and let the volume flow rate in the 
(q, r) channel be V,. If the concentration of 
at the inlet to all channels is varying 
sinusoidally, then, at the end of the (q, r) channel, 
the concentration at any time (varying about 


Channels and 


of different lengths 


There are no 


and 
solute 


the time-average value), will be 


X,, = Ay, 


qr 


cos (wt 4). 


If the liquid issuing from all the channels at 
any time ¢ is rapidly mixed, the concentration 
of the mixture will be given by 


X, = A, cos (wt — ¢,) 


ge Ag 008 (ot — Fer) (20) 
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where L denotes the measured value at the exit 
of a bed of length L, and 


r @ 
V = LL n,, V,, = total flow rate. 


By expanding (20) trigonometrically we have 


and A, sin ¢; = Ay sin ber 


The magnitudes of 4,, and ¢,, for given con- 
ditions will depend upon whether diffusion in 
the channel plays an important part or not. 
If diffusion plays only a small part then, for 
laminar flow in a circular tube, it can be shown 
that 


A, COS dg = 
A, [cos Nr 


Ag, Sin Aj, 


A,, Sin dy, = Ag {sin Ay 


where A, = wl,/[(U)ar], maximum 


velocity (along axis of channel), 


Ci(z) is the cosine integral of 2, 


Si(z) is the sine integral of 2, 


for which values are given in tables of functions. 
These two equations, combined to eliminate 4, 
give the solution to equation (24) of Kramers 
and ALBERDA [6] who give a curve (their Fig. 3), 
showing the variation of A/A, with 2wl U,; 
(present author’s symbols). 

On the other hand, if longitudinal diffusion is 
of importance, then 4,, and ¢,, may be obtained 
by putting £,=1,=—0 (12) 
(remembering that, when /= 0, Q/2w'l=0 
and R/2w'l = 1), to get 


into equation 


sin 


Por 
Ayexp [4/ F’—y(S'/A+ FF" 2)]l, 
/(S'/A— F’/2)]l, 


cos 


sin [ 


(25) 
cos 
in. 2 


= —_ 26 
4(Dg,)* { ) 


(27) 


S' = [FF + (6). 
(cf. The definitions of F*, G* and S* given by 
equations (13), (14) and (15).) 

The longitudinal diffusivity is written as D,, 
because it depends not only on the molecular 
diffusivity but also upon the diameter and length 
of the pore and upon the velocity of the fluid. 
The magnitude of the effective longitudinal 
diffusivity has been discussed by Tay or [8, 9], 
Boswortu [10] and Van Deemrer et al. [11] 
for laminar flow, and by TayLor [12] and Keyes 
[7] for turbulent flow. Levenspre. and Siri 
[13] give a general treatment and McHenry and 
Wituecm [5] deal with a packed bed rather than 
a single tube. 


(28) 


For the sake of developing the argument we 
Using 
the one proposed by Tay.tor [9] (and remember- 
ing that certain conditions must obtain for it to 
be valid), then 


need an analytical expression for D,,. 


p, - Bi Wade 


29 
al 192D (29) 


If we make the 
reasonable assumptions that the pressure drop 


for laminar flow in a circular tube. 
between the ends of every pore is the same and 
is given by Poiseuille’s equation then 
‘ 
Vy 7 R,* p 
qr ? 
l, mé 


where p is the pressure-drop, € is a constant 


(30) 


which depends upon the cross-sectional shape of 
the pores, and pg is the viscosity. Further 
R* He " 
Ny i. ms ’ (31) 


where «¢,, is the fractional voidage contributed 
by n,, channels of length /, and radius R,, and 
H is the total volume of the packed bed. Also 


Nor = po - g, (2k? p)- 


For the convenience of the algebra let 


' 48D & 7 
N, = ( —— F) w; 


d p* 


(32) 
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mM —2*&Pe 
p 


and it will be seen that numerical values can be 
put on M and N. Using equations (21) to (34) 
the amplitude and phase angle at the outlet 
end of the system can be given by expressions 
which involve (a) quantities which are known 
and (b) the minimum number of unknown 
variables. We get, for the case where longitu- 
dinal diffusion is small, 


(34) 


VEA,cos ¢,/(pH Ay) 
[cos A,, — A, 


: x Se, R 


¢ /',') 


sin A, + at Ci (Agr) | (35) 


(similarly for A, sin 4,) and, for the case where 
diffusion is appreciable, 


+ o/ [1 + (NZ/M) B,*))}) 


a: It [1+(N? M*) R,')- 


and these equations can be used to analyse the 
distribution of pore lengths and diameters in a 


a) » (36) 


way similar to that described for the first model for 
the analysis of the distribution of pocket sizes. 
For if w; = L.H.S. of equation (35) or (36), 
and »,,; =the R.H.S. of equation (35) or 
(36), with ¢,, omitted, then, at any value of w,, 


rq 
a y ’ 
Bj x2 €qr"grj 


So certain values of R, are chosen (the order of 
magnitude being known), and for each value of 
R, certain values of /, and for each 
combination of R and | one calculates the value 
of »,,; for a particular value of w; j; contains 
only known or measurable quantities, so the 
simultaneous equations 


are chosen, 


By = 41 M11 7 + €in Mims ete. 


can be written down and solved for «,,. 
As a check, the condition 


rq 
ZS, 


must be satisfied; (« can easily be measured). 


There is also a check on the value of £ chosen for 


r@ 
V = Hp/gp 22 (R,/1,) eg, 
and V can be measured. 


Note that equation (35) will not differentiate 
between pores having the same ratio of length 
to diameter (since 1/R is involved), but equation 
(36) will do so. 

In view of the disadvantage coupled with the 
use of equation (35) as given and because of the 
amount of computation involved by using a 
range of values of R,, it is convenient to use a 
single mean value of R, which seems a justifiable 
step in most cases. 

Then (35) becomes 


> €, (7) [cos A, — A, sin A, + A,? Ci(A,)] (87) 


for the case where diffusion is small, 
A, = péw 1? /(2pR*) and (86) becomes 


itp I! = 08 $y E: (i; i) 

exp [ Ml, 1+ 4/[1 + (N2/M*)R*) 
5 Fe eee 

TR a (a/ [* + - -1)}] (9 


for the case where diffusion is appreciable. 
Then we have 


where 


r 
Bj — xe, Nrj 


which, on writing out in the form of simultaneous 
equations for various values of w,, needs fewer 
terms in the sequence. 


The checks will be as before. 


Some FurtTHER CONSIDERATIONS 


In the treatment of model 1 it was necessary 
to put a value on U. This may be done either by 
(a) arranging the flow in the channels to be 
laminar and such that diffusion plays only a 
small part. A value of R/L, for the (identical) 
channels may be determined by using equation 
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(87). If a value is put on L, (by assuming a 
“tortuosity factor”) and hence on R, then, 
since 


V = Cp w R*/(L, né), 


where V = volume flow rate, 
€ = constant; usually between 2-5 and 
3-0 [1], 
C = number of channels, 


L, = 


. = effective length of the channel, 


a value can be put on C and hence on U, 

Since the velocity at the walls is zero, the 

concentration does not change at the walls or 
beyond; i.e. the pockets do not contribute. 
This method assumes that the fractional volume 
of dead space is unaltered when the flow in the 
channels alters. 
(b) making o large; for, from equations (7), (8) 
and (9) it will be seen that if o is made large, 
Q/o and R/o become small and the pockets 
again will not The system then 
approximates to that of model 2. In this method 
the velocity in the channels is unchanged. 


contribute. 


CONCLUSION 


A method has been suggested whereby a flow 
structure or a physical structure may be analysed 
in terms of models which have to be chosen so 
that they approximate to the system under 
investigation and whose effect on a flow of 
fluid containing a sinusoidally varying concentra- 
tion of solute may be mathematically analysed. 
The attenuation and phase angle of the output 
streams (assumed rapidly mixed in a small mixer) 
are measured at a number of different frequencies 
and for each frequency an equation, derived from 
a consideration of the effect of the model on the 
sine-wave, is set up in the form given, for example, 
by equations (19.1) to 19.n); one of the co- 
efficients (7,;,) contains r-dependent terms, the 
values of which are known (or assumed and 
afterwards checked) and j-dependent terms. 
The unknowns (8,) are what it is wanted to find. 
The constant for the equation (4;) contains only 
terms which are known or measurable at the 
particular frequency, «;. 

The number of values of the frequency is 
equal to the number of values that are given to 


the r-dependent part of »,;, so that the 
simultaneous equations are soluble. 

The decision as to the number of values to 
give to r depends upon the closeness of the 
approximation desired, the upper limit being 
determined by the amount of computation or 
by experimental limitations. 





" y 
S28! 


(7-- 


Fic. 5. 


Typical histograms of the distribution of relative 


volumes 8, of pocket length J,. 


As an illustrative example, Fig. 5 shows what 
might be the actual distribution curve (volume 
against depth) for pockets associated with a 
channel, and the histograms obtained by making 
three at these 
numbers of frequencies, and solving one, three 


one, and ten determinations 


and ten simultaneous equations respectively. 


Results have been given for models which 
illustrate the method and which are likely to be 
useful in practice, but the choice of model and 
the interpretation of the results will depend 
upon the problem in hand. For example, it 
might be suspected that a catalytic reactor has 
an appreciable fraction of its void space “* dead ” 
i.e. filled with stagnant fluid. An analysis based 
on model 1 would give the fraction of dead-space 
as if the reactor were actually of this model form. 
In practice this will not be true; e.g. the dead- 
spaces will not be of prismatic cross-section, but 
the analysis will give some useful quantitative 
results, which would be even more useful if they 
were applied to the comparison of two roughly 
similar reactors. 

On the other hand the reactor may be suspected 
not of having dead spaces but of having channels 
in which some of the fluid spent too long, and 
others in which it spent too short, a time. An 
investigation, based on model 2 would give a 
distribution of channel lengths and diameters as 
if the reactor were of this form. 


As pointed out above, the model can be 
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changed, even for the same reactor (for example), . = Length of channel (model 1) 
by altering the conditions; viz. frequency or ve = Effective length of a channel 
velocity of fluid, but care must be taken that the Hydraulic radius 
system does not change also. ‘ 246 Du/p 

The choice of magnitudes (in particular of ij = 48 DE? wo; ?/p? 
velocity), will depend upon the physical quantities, (: ah Ze'y + cos Se "’ 
which would in general be known roughly (e.g. cosh 2a'f + cos 2a | 
channel diameter), and upon the conditions 
desired. For example it might be desired to use 
an analysis in which diffusion in the channels ae 
was important, and a velocity would have to be awe + eaSws 
chosen to make this so. 

In general it would be true that the nearer the 
chosen model approached actuality, the more 
difficult the mathematical analysis and the more 
laborious the calculations. 


Pressure drop 


Average concentration in a pocket 


sinh 2a 7 + sin 2w'l 
- — (model 1) 
cosh 2a'l cos 2a 
Radius of a channel (model 2) 
v [FP + @] 
y (CPt? + (G*)?] 
VPP +(e?) 
Sine integral 
Acknowledgements—The statement of the problem Time 
and suggestions for the theoretical and practical ' Average velocity 
methods of attack are due to Dr. J. A. Srorrow. ', = Maximum velocity 
The author is indebted to his colleague, Volume rate of flow 
Mr. O. E. Porrer, for very helpful suggestions N ~ Concentration relative to a mean value 
and criticisms. Distance from closed end of a pocket 
Distance along channel 
Volume of pockets in unit length of channel 
NOTATION n/a 


= Cross-sectional area of a channel Voidage 
= Amplitude Q/2 w 1 (model 1) 
Complex amplitude = Ae? » = Defined by equations (35) or (36) et 
Number of identical channels (model 2) 
Cosine integral Phase angle 
Generalized longitudinal diffusivity ‘ wl /2U 
Molecular diffusivity 2S £,Q,/(2«'l,) (model 1) 
; VéAy (pH Ag) cos OL (model 2) 


‘u\* — Bo 
(<) + 2Dw'l Dynamic viscosity 
Constant 


U . Ls 2w 1 
») : Phase angle 


(U 2D/* Angular velocity 


\ (@/[2D)). Dimensions : L 


1— BR [2w'l)) 


7 


n Identifies a value of w 
(: - 2 8, R, (2, «) L Value at exit of bed 
1 j 


Value at entrance to bed 


@w 
D 


= Identifies a radius of a channel (model 2) 


D Identifies length of a pocket (model 1) 
H = Total volume of bed Identifies a length of a channel (model 2) 
| = Length of a pocket (model 1) ' Value at distance y from closed end of pocket 


| = Length of a pore (model 2) = Value at distance z along channel 





The flow-structure in packed beds 
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Correlation of scraped-film heat transfer in the votator 
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Research & Development Laboratories, Thomas Hedley & Co., Limited, Newcastle-upon-Tyne 
( Received 20 May 1957) 

Abstract This work on a votator was undertaken to correlate the scraped film coefficient of 

heat transfer, /,, with operating conditions, plant dimensions and physical properties of the fluid. 


All the available data, which include runs on the cooling of glycerol, water and two similar 
glveeride oils in four votator-assemblies, are provisionally correlated by the dimensionless 


hed, 1. Dg ‘" Cp Hs - 0,17 men 
\, .; h v L 


The ratio of shaft to tube diameter, D,. D,, varied little throughout the work, from 0-75 to 


equation : 


0-81, so that any effect of this ratio is incladed in the constant of 4-9, 

There is some indication however, of the existence of flow regimes analogous to streamline, 
transitional and turbulent conditions in flow through conduits and in conventional liquid mixing 
operations. This evidence indicates that exponents on the groups D,N/o and Dp p/m are 
higher in the turbulent range than the values appearing in the above equation. 


Résumé Ce travail sur un “ votator,” a été entrepris dans le but de relier le coefficient 
d’échange de chaleur a travers le film balayé, 4,, aux conditions opératoires, propriétés physiques 
du fluide et caractéristiques de lappareil. Toutes les donnés, disponibles, comprenant des essais 
le refroidissement du glycerol, de Peau et de deux glycerides homologues dans quarte assemblages 

de votators, satisfont a la relation sans dimension suivante : 

0,57 0,47 )\ O17 
AD, 9 Dp °) (° “) DN ”) 
hk m hk v L 

Le rapport de arbre au diamétre du tube D,/D, a peu varié au cours de ces travaux (de 0,75 
& 0,81), de telle sorte que la constante 4,9 suflit & rendre compte du phénoméne pour cette faible 


0,37 


variation. 

Il existe cependant quelques indications sur les régimes d’écoulement analogues aux régimes 
laminaires transitoires et turbulents dans les conduites et dans les opérations de mélange. 

Ces remarques montrent que les exposants des groupes D,N/o et Dw p/ sont plus grands en 
régime turbulent que les valeurs mentionnées dans l’équation précitée. 


Zusammenfassung—-Diese Untersuchung an cinem Votator wurde unternommen, um die 
Wirmeiibergangszahl beim Abstreifen, 4,, mit den Arbeitsbedingungen, den Abmessungen der 
Anlage unde den physikalishchen Ligenschaften der Flissigkeit in Bezichung zu setzen. Alle 
verfiigbaren Daten, in denen Versuche tiber die Kihlung mit Glyzerin, Wasser und zwei dholichen 
Glyzerin-Olen in vier Votatoranordnungen enthalten sind, kOnnen vorliufig durch folgende dimen- 
sionslose Bezichung dargestellt werden. 

7) 

L 


h,D, as (7° ‘ (‘n)"" (2) 
k » h v 

Das Verhiltnis von Schaft —- zu Rohrdurchmesser (D,/D,) dindert sich nur wenig innerhalb 
dieser Untersuchungen, nimlich von 0,75 bis 0,81, so dass ein etwaiger Einfluss dieses Quotienten 
in der Konstaten 4,9 enthalten ist. 

Es gibt einige Anzeichen fiir die Existenz von StrOmungsbereichen Ahniich der laminaren, 
Obergangs - und turbulenten Form bei der Kanalstromung oder den iiblichen Mischverfahren 
von Flissingkeiten. Im turbulenten Bereich sind niimlich die Exponenten der Kenngréssen 
D.N/o und Dp p/p hoher als sie in der obigen Gleichung erscheinen. 


0,17 0,37 
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Correlation of scraped-film heat transfer in the votator 


INTRODUCTION 
Tue Votator, a double-pipe, scraped surface 
heat exchanger, is widely used in the food and 
chemical industries, for example in the erystalliz- 


ing of petroleum waxes, the plasticizing of mar- 


garine and shortenings and in sulphation and 
sulphonation reactions. 

Despite its widespread applications, only one 
paper has been found dealing quantitatively with 
heat transfer in the votator. This paper, by 
Hovurron [6], deals with a water-to-water heat 
exchange in a small laboratory unit, in which 
shaft speed, jacket water velocity and throughput 
rates were studied in relation to the over-all and 
film heat transfer coeflicients., 

No correlation of the scraped film coetlicient of 
heat transfer with operating conditions and 
physical properties has previously appeared, 
although the need for such a correlation was 
underlined thirteen years ago by Hovu.ton, 
whose data is used in the present work. 

Many of the 
incidentally during work on a more extensive 
investigation. If the study of heat transfer 
through the scraped film had been the prime 
object of that work, cooling would have been 
effected by brine rather than by direct-contact 
ammonia evaporation, thus simplifying measure- 
In an ideal case 


present data were obtained 


ment of the heat transferred. 
also the fluids on the scraped surface side would 
have been chosen to avoid phase change. 


EXPERIMENTAL WorK 


-Ipparatus 

Two ammonia-cooled votators were used, they were 
structually similar and the small unit is shown in Fig. 1. 
They were designed by the Girdler ¢ orporation (ULS.A.) 
and made by A, Johnson & Co, (London) Ltd., the details 
being as follows 


Table 1 





Small 
V otator laboratory, Large unit 


unit | 








| ors 
| 7-99 11-08 
| #875 4375 


No. of tubes in series 
“ Effective ” 
Heat transfer tube o.d. (in.) 
Heat transfer tube i.d. (in.) | : nO 4-00 
Revolving shaft o.d. (in.) $25 3-25 
A. Johnson & Co, specification no. | 3/414 M1234 


scraped area, A,(ft*) | 


2-25 





Attached to each shaft and 180° apart were two scraper 
blades of the “ floating * type, which moved freely to the 
wall under centrifugal force and fluid kinetic pressure as 
the shaft revolved. (The shaft plus blades is called the 
mutator). 

The votator tubes were made of mild steel with thinly 
chromium plated inner surfaces. 


Materials used and their physical properties 


Water and glycerol, Since the glycerol contained a small 
but unknown quantity of water its viscosity was measured 
in a Ferranti viscometer ; the results are shown in Fig. 8. 
bor the other physical properties those for pure glycerol 
were used. Two similar glyceride oils, A and B, were also 
used and the following typical values were taken for their 
physical properties : 


Specific heat: (partially crystalline oil beyond 
votator) 0-45 B.t.u. /lb °F 


Specific heat : (less erystalline oil within votator) 


OS B.tu. ‘ib °F [3] 
Latent 95 B.t.u. Tb. 


0-1 B.t.u. hr ft? -F/ft [3] 


heat of fusion : 
Thermal conductivity : 


Specilic gravity : O-sS®%. 

Progressive crystallizing of the oils along the votator, 
plus substantial super-cooling, prevented accurate assess- 
ment of an average viscosity within the device. 

Values in Fig. 9 for similar oils at temperatures above 
cloud point [9] were extrapolated through about 20°F to 
the 
\ reason for accepting this procedure is given later. 


provide values within relevant temperature range. 


at a given 
estimated * Solid 
(S.C.1.) using the method and termin- 


The equilibrium percentage solid content 


temperature experimentally as 
Content 


ology of Fuuron, Lerron and Wiixe [4]. 


was 


Index ~ 


Measurements were as follows : 


Table 2 





Ou A 


14-5 





Linear interpolation has been found to give a good 


approximation between these temperatures for these 


oils [9]. 


Determination of the experimental scraped film coefficients 


Experimental conditions and results from runs on the 
small votator appear in Table 3, values calculated from 
HovuLron’s runs are given in Table 4, while values ob- 
tained on the large votator are given in Table 5. 
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insulation 


Section ot A-A 


Longitudinal and cross-sectional diagrams of 
ammonia-jacketed laboratory votator. 


Over-all coefficients based on the inside scraped surface 
area of the votator were first caleulated from the general 
equation : 

Q = UA, LM.T.D. 


effective - was calculated 
in accordance with Houron’s assumption [6] that the 
flanged heads of the refrigerant half 


effective for cooling. 


where the “ seraped area, -1,, 


jacket were only 


In the case of the glyceride oils A and B, the heat 


transferred per hour Q, consists of : 
Heat equivalent to about 85 per cent of the energy 
consumed by the mutator motor | 
Sensible heat due to cooling the partially erystall- 
ized oil. 
Latent heat resulting from some crystallization of 
the oil in the votator. 


The fraction crystallizing in the device was calculated 
as the difference the S.C.1. at the 
temperature reached by the partially solidifying oil after 


between maximum 
emerging from the votator and the crystallization cor- 
responding to the temperature rise beyond the votator. 
Allowance was made for frictional and mechanical heat 
effects beyond the votator from preliminary heat balance 
measurements. 

Attempts to calculate the ammonia film coeflicient 
either by assuming a forced convection mechanism and 
using the usual Dittus—Boelter equation, or by assuming a 

usual 
Prandtl 


natural convection mechanism and using the 
correlation involving Nusselt, 


numbers, led to an ammonia film coefficient less than the 


Grashof and 


overall coeflicient—-an impossible result. 

It was thus concluded that a much higher ammonia 
film coefficient resulting from a nucleate boiling mechanism 
must be prevailing. This was partly confirmed by visual 
observations of vapour plus entrained liquid ammonia 
flowing upwards past a glass window in the upper outlet 
pipe from the ammonia jacket. 


Ammonia film coefficients were therefore calculated 
from a correlation by McNexry [7], based on a nucleate 
boiling mechanism. The properties of the ammonia film 
were evaluated at the mean film temperature. If either 
the bulk temperature or the adjacent metal wall tem- 
perature had been used the 
would have differed from that corresponding to the mean 


resulting film coefficients 


film temperature by less than 2 per cent, which is well 
within the limits of accuracy of McNewty’s correlation. 
Seraped film coefficients were than extracted from the 
overall coeflicients by the general equation : 


l 1 1 
— = +—+ 
U hyn, he 


1 
h, 


Here h,, is the usual metal wall coefficient, given by 
dividing the thermal conductivity by the metal thickness. 


DIMENSIONAL ANALYSIS 


Heat transfer through the scraped film will be 
by conduction the thermal con- 
ductivity, k, of the fluid is involved. The film 
thickness should depend on the linear axial 
velocity, v, of the fluid, on fluid density, p, 
viscosity, , tube diameter, D, mutator shaft 
speed, N, shaft diameter, D,, and upon the 
diameter of the mutator, D,,, which last is defined 
as twice the distance from the blade edge to the 


and hence 


shaft axis. In view of the phase-change of some of 
the oil, and the possibility of modifications in 
flow pattern along the freezer, the tube length, 1, 
and blade length J,, may be influential factors. 
Also, for a given rate of heat transfer, the bulk 
temperature of the stream is affected by the 
specific heat of the fluid, c,, which must also be 
included. 

The dimensions chosen as fundamental in the 
length, time, heat and 
temperature, necessitating the inclusion of the 
mechanical equivalent of heat, J — a dimensional 
constant. Thus the general equation relating h 


analysis were mass, 


s 
to these twelve factors in exponential form and 


where ¢ is a constant is : 

h, = (Dg D? Df Uf Nf v8 ef kp" uw? J®) (1) 
Dimensional analysis of this relationship leads 

to the general equation : 


b 


Poe CePA GY 
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l, almost equals l,, equals L (say), r equals d + e, 
q was found to be zero and the assumption has been 
made that D,,/D, is unity, a point commented on 
later. 

The Nusselt and Prandtl groups have their 
usual significance, while the group Dy p/u only 
resembles the Reynolds number since it contains 
the tube diameter, D,, instead of the equivalent 
diameter of the annulus between tube and mutator 
shaft. The group D,N//v is the ratio of a measure 
of the rotary peripheral velocity to axial velocity 
of the fluid and the groups D,/D, and L/D, 
partly define the geometry of the votator. 
Powers of the dimensionless groups and the 
coefficient 

Physical properties have been evaluated at 
average bulk fluid temperatures : 

Exponent f was evaluated by using the data for 
oil A and glycerol in log — log plots of h,D,/k vs. 
DN /e. 

Here Dv p/p remained constant for each curve, 
c, #/k remained constant for each substance, and 
D,/D, and L/D, remained constant throughout. 

The slope of Fig. 2 shows f to be about 0-17. 
The significance of the line for water is discussed 
later. 

10000, 


Fic. * 


Exponent n. In a preliminary assessment of 
this exponent, data for glycerol and oil A were 
used in log — log plots of 


Dp p 


AD (DIY oy 
p 


k \o 
where c, 4/k remained constant for each curve, 
and D,/D, and L/D, remained constant through- 
out. 

The slope of Fig. 3 shows n to be about 0-52 ; 
this value was slightly revized later from a final 
plot of all the available data. 








> = Glycerol 
=O1WA 

















Exponent r. Data for oil A on the small 
votator were used with data for oil B on the large 
votator when either two or three tubes in series 
were used. This provided a range of tube lengths 
from 0-87 to 11-44 ft. Log —log plots were made of 


O17 052 L 


(72*) 5 


k rm 


assuming the same c, u/k for these two similar 
oils. 


h,D, (Px ) 


v 
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The ratio D,/D, was no longer strictly constant, 
but the variation is considered small, i.e. for the 
small votator, D,/D, was 0-75, while for the 
large unit it was 0-81. 

The slope of Fig. 5 shows r to be about — 0-37. 

Exponent i. For glycerol, oils A and B and 
water, log — log plots were made of 

AD, (De)*™ (DiN)*™ (D)*" 5 ot 

k mm 2 L k 

A value of 0-47 was provisionally taken for i 
from the slope of the line in Fig. 6. 

Exponent b. The effects of varying the ratio 
D,/D, were not studied and it remained nearly 
constant throughout both this work and that of 
Hovutton. The factor (D,/D,)” is therefore 
included in the value of the coefficient found 
below. 


Coefficient ¢ and revised value for exponent n. 
Because of the wide range of values of Dw p/ 
covered in this work a final plot of all the data 
against this group was made to assess the co- 
efficient ¢. Fig. 7 shows a log — log plot of : 


MD (Pe yon (een) (7 ve, Dee 
k v, ( k ) L in 
The intercept has a value of 4-9 and the slope of 
the line gives a revized value for n of 0-57. 

1000 





Slope =c.a. -O-37 


| 
O Oil A, smal! votator 
© Oil B, large votator (2 tubes) 
+ Oi1 B, large votator (3tubes) 


-0 "exw? 


mn 


i) 


*'C 
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Slope = c.a. 0-47 














C Glycerol. small votator 

© Oil A. small votator 
Reemenduia © Oil B. large votator. (2 tubes) 

+ Oil 8B. large votator. (3tubes) 

% Water. smal! votator 

x Water. (from Houlton) 


| 
| 
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The equation relating to the line on Fig. 7, 
which contains all the available data, is :— 

hD, 

k 
Dw p 0-57 Cy ft 0-47 DN O17 D 0-37 

9 t p t t B 

so (~—) () () G&G) © 


v 


The equation has been determined for values of , 


D,/ D, between 0-75 and 0-81, 


the 
turbulent flow regimes 


Evidence for existence of streamline and 

It is well known that the transition between 
streamline and turbulent fluid flow through 
straight conduits is characterized by values of 
Dw p/p between 2000 and 4000, while the 
corresponding transition region for conventional 
liquid mixing vessels is defined by values of 
D2ZN p/p between 10 and 10,000 [5]. 

Maximum and minimum values of these two 
forms of Reynolds number are shown for these 
runs in Table 6, where the equivalent diameter 
of the annulus between tube and mutator shaft, 
D,, is given by D, — D, and where the impeller 
diameter, D,, is taken as D,,. From this table 


analogy with conventional liquid mixing vessels 


1000 


would indicate that turbulence is only fully 
developed in the case of the water runs. 


Table 6. Maximum and minimum values of 


(D, — D,)v p/p and DZNp/u 





Fluid (D, D.)v p/p D,2N p/p 


meni * 
Min. Min. 


Maz. Maz. 


Glycerol 
Oil A 
Oi B 
Water 


0-7 1-4 
37 O-4 
63-0 135-0 
505-0 1690-0 


MeO | 
134-0 

1690-0 | 
28200-0 


175-0 
670-0 
2030-0 
380,000-0 





There is also some indication of the existence 
of two regimes of flow in Fig. 2, where the middle 
four values for water from Table 4 are plotted to 
indicate an exponent of the order of 0-5 for the 
group DN vin the region. This 
exponent cannot be with greater 
the available data (that of 
Hovutron) since, although the values of L/D, 
D, D, and ¢, wk corresponding to points on this 
line remained 


turbulent 
determined 


precision from 


constant, values of 


ranged from 5910 to 6350, 


D,v p/p 





Slope =0-57 
O=4-9 


0-47 < 0-37 
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Fig. 4 suggests an exponent of the order of 0-9 
for the group Dw p; in the turbulent region. 
Although the groups L./D,, D,/D, and c, w/k are 
constant for Fig. 4 more data are clearly needed 
to establish this exponent under turbulent 
conditions. Analogy with correlations for stream- 
line and turbulent flow through conduits indicates 
that the exponent on the group LD, will at least 
be reduced for turbulent due to 
reduced entrance effects. 

It seems from the above observations that the 
exponent on the Prandtl group in equation (3) 
may be a hybrid value obtained by combining 
data for streamline and turbulent flow. Accord- 
ingly the correlation regarded as 
provisional pending further data on the existence 
of more than one flow regime and on the various 


conditions, 


must be 


group-exponents at high Reynolds numbers. 

It appears, then, that the situation regarding 
equation (3) may resemble that underlined by 
comments of Prarr [8] when comparing his 
own correlation that of Cumrox, Drew 
and JrBENS [2] for heat transfer in agitated 
vessels. Pratt's comments strongly suggest that 
the well known single-line correlation of CuiLron, 
Drew and Jesens actually spans more than one 
Comments relevant to this latter 


with 


flow regime. 


topic also appear in discussions of Brown, 
Scorr and Toyne [1]. 
Features of equation (3). The constants in the 


correlation have been found from experiments in 
which the fluid was being cooled, The fluid adjacent 
to the scraped metal wall would thus be colder 
and more than in the 
Slightly different constants might be expected 
if the fluid were being heated, since the tem- 


viscous main steam. 


perature gradient and hence the viscosity pattern 
reversed, 

A comparison between the exponents on the 
groups D.N/v and Dw» p/ indicates that the 
rate of heat transfer is not dominated by the 
mutator speed and that the contribution afforded 
by linear flow velocity through the votator must 
be considerable. 


would be 


Logarithmic mean temperature difference 


In calculating over-all coefficients from the 
experimental data the logarithmic mean tem- 


perature difference between ammonia and the 
votated fluid was used ; this practice is only valid 
if constant specific heat can be assumed as in the 
case of water and glycerol. It is believed that this 
assumption can be largely justified in the case of 
oils A and B despite their partial phase change 
within the votator, in view of their substantially 
linear relationship between equilibrium solid 
content and temperature. Thus, over most of 
the temperature range through which the oil is 
cooled, it has an effective specific heat com- 
pounded of nearly constant amounts of latent 
and sensible heat. 


Physical properties of oils A and B 

A perfect correlation in Fig. 7 is partly prevented 
by the shortage of accurate physical data for 
oils A and B, necessitating extrapolations and 
the use of typical values. 

In particular, since the true viscosities for oils 
A and B at sub-cloud-point temperatures are 
not known, it may be felt that the data for these 
substances lie on the general line of Fig. 7 by a 
coincidence connected with the use of “ pseudo "’ 
viscosities derived from the extrapolation of 
Fig. 9. However, it should be noted that ex- 
ponents on groups containing viscosity (namely 
Dp p/p and c, uk) have not been derived ex- 
clusively from data for oils A and B, but have 
been obtained independently from data for 
glycerol in the case of Dy p/p (Fig. 3) and from 
data for glycerol and water in the case of c, 4/k 
(Fig. 6); true viscosities for glycerol and water 
are, of course, known. 

Furthermore, partial cancellation in the groups 
Dy p/p and c, u/k leaves viscosity to the rather 
low power of — 0-1 which may partially account 
for the close agreement of data for oils A and B 
with the provisional correlation, despite the use 
of “ pseudo” viscosities in place of true values 
for these substances. 


Contact between blades and tube wall 


It is known that even slight clearance between 
the edges of the scraper blades and the tube wall 
may lead to considerable reductions in heat 
transfer [9]. Because of the “ floating” blades 
used in this work the assumption has been made 
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that close contact between blade-edges and wall 
is maintained; this should be true at least for 
water and glycerol, and by taking D,,/D, equal 
to unity throughout the assumption has been 
extended to the partially crystallizing oils A 
and B. 


Factors not included in the correlation 


The effects of variation in mutator design have 
not been included, for instance with regard to 
the shape of the cross-section of the shaft, the 
variation of the ratio D,/ D,, the number of blades, 
and the clearance between blade-edges and tube 
wall. More work is needed to incorporate the 
above factors and to distinguish more clearly 
between streamline and turbulent conditions. 
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NOMENCLATURE 


“Effective” scraped heat transfer surface ft? 
Specific heat B.t.u/ib°F 
= Conduit inner diameter, or equivalent 
diameter of an annulus ft 
Impeller diameter 
Mutator diameter (2 x distance from 
blade edge to shaft axis) 
Diameter of mutator shaft 
Internal diameter of Votator tube 
Metal wall of heat transfer 
B.t.u. /ft®hr°F 
= Ammonia film coefficient of heat transfer 
B.t.u. /ft®hr°’F 
= Seraped film coeflicient of heat transfer 
B.t.u. /ft@hr°F 
Thermal conductivity B.t.u./ft®hr°F /ft 
Length of scraper blade (= L) ft 
l Length of votator tube (= L) ft 
L.M.T.D. 


coefficient 


Log, mean temperature difference 
between fluid being cooled and refrigerant °F 
\ Mutator shaft speed rev /hr 
Q = Total heat transferred/unit time B.t.u./hr 
u Over-all coeflicient of heat transfer 

B.t.u. /ft®hr’F 
Average axial flow velocity of fluid 
through votator 
Fluid density 


ft/hr 
ib /tt® 
Absolute fluid viscosity Ib /ft hr 


A proportionality constant 


Exponents : 
a, b. . d, ¢, f. SS & J," P. q- 


r d +e. 
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REPORT OF MEETING 


Symposium on Chemical Reaction Engineering 


Syllabus 


Chemical reaction engineering is a part of chemical engineering in general. 


It is a 


new branch of science which is still in the development stage. It aims at controlling the chemical 


conversions on a technical scale and will ultimately lead to appropriate and successful reactor 


design. 
transfer, and reaction kinetics 


An important part is played by various factors such as flow phenomena, mass and heat 
It will be clear that in the first place it is necessary to know 


these factors separately. Yet this knowledge in itself is insufficient. The development of chemical 


conversions on a technical scale can only be understood from the relation and interaction between 


the above-mentioned factors. 


On 7, 
posium on Chemical Reaction Engineering took 
the Royal 
Institute. This symposium, which at the same 


8 and 9 May 1957 First European Sym- 


place in Amsterdam in Tropical 
time was the twelfth conference of the European 


Federation for Chemical Engineering, 
ganized by the Chemical Engineers Section of the 
Koninklijk Instituut 
Chemical Technology and Industrial Chemistry 
Section of the Koninklijke Nederlandse Chemische 


Vereniging. 


was ofr- 


van Ingenieurs and the 


The members of the Symposium Committee 
were: Dr. Ir. J. C. Vuiwcrer (Chairman), 
Prof. Ir. E. F. Boox, Dr. C. van Herrpen, 
Dr. Ir. A. Kurnkennerc, Dr. K. Rierema 
(Secretary), and Dr. Ir. G. S. per ViuiEs 
(Treasurer). 

There were three official languages : 
Three 


VAN 


French, 
German and English. the 
symposium the members had received preprints 
of the fourteen papers 
were dealt with in five successive sessions. Each 
session had a chairman as well as a rapporteur, 


weeks before 


contributions. These 


whose task consisted in reporting on the discus- 
sions and the session. 

In all, 260 members had registered for the 
symposium, half this number being Dutch, the 
other half from twelve different 
The fact that there were three languages has not 


countries. 


caused any difficulties, which is also due to the 
good oflices of the linguistic adviser, Mr. W. G. 
pE Borer. 

The symposium was opened by Dr. ViucrTER, 


This relation and interaction will be the main theme of the 
symposium, “ Chemical Reaction Engineering”. 


who welcomed the official guests and drew atten- 
tion to the pioneering work done on behalf of the 
European Federation for Chemical Engineering 
by Mr. Jean Géranp, who recently died in Paris. 
Dr. Viccrer added a special word of welcome 
to the students, twenty-five of whom attended 
the symposium. Finally he expressed his conti- 
dence that the symposium would be a success. 


First Session: “ INrropuctrory Papers ” 


Chairman Dr. Ir. J. C. Viwerer. 


Rapporteur Dr. Ir. D. va~N Zoonen. 


Speakers Prof. Dr. D. W. van Kreveten (Central 
Laboratorium, Staatsmijnen in Limburg) ; 
Prof. M. Lerorr (Faculté des Sciences de 
l'Université de Nancy) ; Prof. Ir. H. Kramers 
(Laboratorium voor Fysische 


Technische Hogeschool, Delft). 


Technologie, 


In this session introductions to the symposium 
were given. The first was a general introduction 
by Prof. van Kreve.en, who explained the aim 
of the symposium and gave a detailed classifica- 
tion of the field of study. 
shown how the various papers to be presented at 
the symposium fitted in this classification. 

In a historical survey Prof. van Kreve cen 
showed that three periods can be distinguished 


In addition, it was 


in the development of this branch of science. 
The first is the “ empirical-phenomenological ” 
period up to about 1915; then 
* molecular-statistical "’ and, from 
“ technological ” period. 


comes the 
1937, the 
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One of the most important aspects of the 
technology of chemical kinetics. 
This divided into microkinetics 
macrokinetics. Microkinetics deals with chemical 
kinetics or molecular-scale processes. To these 
belong : spontaneous reactions, catalysed reac- 
tions, chain reactions, as well as side reactions 
and Often it is 
necessary, for the design and calculation of a 
chemical reactor, to know the chemical kinetics 


reactors is 


can be and 


consecutive reactions. not 


of all the stages of the reaction mechanism. 
In general it suflicies to have a general picture 
of the kinetics. 

In the laboratory it is usually possible to make 
the reactions proceed under ideal conditions and 
to choose dispersion, homogenization as well as 
heat supply and removal at will. In industry 
often The 
complicated — by transport, 


this is impossible. problems are 


then mass heat 
transfer and transfer of impulse on a larger scale 
than the molecular one. This is called macro- 
kinetics, which deals with phenomena proceeding 
that 
the dimensions of the dispersed 


the thickness of the 


on a scale is comparable with reactor 
dimensions, 
substance or 
layers. 

The influence of macro- and micro-kineties is 


boundary 


reflected in the way in which the reaction rate 


depends on temperature and pressure, in the 
stability of the reactor as well as in the yield 


and selectivity of the reactions. 

Prof. VAN KREVELEN the 
words: “I think I may say that in this field of 
study, with its numerous aspects, our symposium 


concluded with 


is the most comprehensive survey ever given”. 

In the second introductory paper Prof. Lerorr 
gave a survey of microkinetics and in particular 
of the theory of the chemical kinetics of homo- 
geneous gas-phase reactions. 

The most important quantity in kinetics is 
the rate. By definition this is the 
quantity of substance converted per unit time. 


reaction 


The reaction rate can be measured by analysing 
samples taken at certain intervals. Preference 
should be given, however, to physical methods 
which enable the course of the reaction to be 
followed by effecting measurements in the reactor. 


Although it is not always necessary to know the 


entire reaction mechanism it is desirable to do 
so for purely scientific work. 
There are simple reactions, 


which may be 


mono- bi- or tri-molecular. As a rule they are 
bimolecular. In addition, there are complex 
reactions, which may be accompanied by side 
reactions or consecutive reactions. Besides, 
chain reactions often occur. They can be classified 
as processes initiating, propagating or terminating 
the chain. These reactions are characterized by 
the fact that one of the links form a_ particle 
identical to one that has disappeared before. 
This was explained by referring to the example 
of the thermal decomposition of acetic aldehyde, 
on which Prof. Lerorr himself has worked for 
several years. It is often very difficult to find the 
mechanism of chain reactions, 

The order of a reaction is the sum of the ex- 
the the 


reaction rate is proportional. Simple mono-, 


ponents of concentrations, to which 


bi- or tri-molecular reactions are first-, second- or 
third-order reactions. If a reaction proceeds 
stepwise or is a chain reaction the order is often 
a fraction. Sometimes the order changes during 
the reaction under the influence of the products 
formed. 
According to Arruentus, the influence of 
temperature on the reaction rate is given by the 
activation This is, so to speak, the 
energy required to lift the reacting molecules 


energy. 


over a potential barrier. In general, the influence 
of temperature on the reaction rate is considerable. 
Finally, Prof. Lerort made a few remarks on 
photochemical reactions. 

Prof. KraAMERS gave an introduction to the 
chemical 


physical factors playing a part in 


reaction engineering. First of all he remarked 
that, for the physicist, it is gratifying to see that 
physical phenomena greatly assist in making a 
Such 
reactants 


reaction feasible. phenomena 
include the 


products to and from the reaction zone, the 


chemical 
transportation of and 
mixing and dispersion, diffusion, heat and mass 


transfer, which all belong to macrokinetics. 
They were elucidated in the paper with the aid 
of some examples worked out by Prof. Kramers 
and his co-workers. 


The first two examples had been taken from 
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the field of hydrodynamics. The first concerned 
the distribution of residence times in tubular 
reactors in the case of laminar flow, turbulent 
flow or flow through a packed bed. The second 
example showed that in a continuous process 
with known residence-time 
impossible to calculate the conversion unless 


distribution it is 


the reaction is a first-order one. For reactions of 
different order the result of the 
depends on the model chosen to explain the 
residence time distribution. For actual practice 
it is a fortunate that in 
cases the differences are not very large. 
In the next two examples, which referred to 
mass transfer, Prof. Kramers showed that if a 
mixture is uniform on the macroscopic scale this 


calculation 


circumstance most 


need not imply that there is homogeneity on the 
molecular scale. It is possible to ascertain 
uniformity, for instance, with the aid of an electric 
conductivity cell and molecular homogeneity 
by following the course of a very fast chemical 
reaction. Another example dealt with 
transfer in a fixed bed. It appears that the data 


from the literature can be made to fit in with one 


mass 


another over a wide range by applying to the 
usual correlations some corrections that are easy 
to justify. 

Finally, two examples were given that referred 
to heat transfer. It was demonstrated that in a 
fixed bed the ratio of grain diameter to tube 
diameter may greatly influence the transfer of 
heat. For a fluidized bed the variables governing 
heat transfer were divided into two 
those referring to thermal quantities and those 
referring to mixing in the bed. The transfer of 
heat to a small object in the bed, to a low vertical 


groups : 


wall and to a high vertical wall must be considered 
separately. 


Seconp Session : * TRANSPORT PHENOMENA IN 


HETEROGENEOUS REACTIONS ” 


Chairman Prof. Dr. D. W. van Kerevecen 


Rapporteur VP. J. Horrizer. 

Prof. Dr. E. Wicke (Institut fiir Physikalische 
Chemie der Universitat Hamburg) ; Dr. Ir. J. 
G. van pe Vusse (Koninklijke /Shell- 
Laboratorium, Amsterdam); Ir. R. A. T. O. 


Speakers 


Nustnc (Laboratorium voor  Fysische 
Technologie der Technische Hogeschool, 
Delft); also on behalf of Prof. Ir. H. Kramers, 


The papers presented at this session all dealt 
with cases in which the velocity of the process 
under consideration was dependent, not only on 
the chemical reaction rate proper, but also on the 
velocity of a mass transport process. In Prof. 
Wicke’s paper this was the case in heterogeneous 
gas reactions ; the other two papers referred to 
the absorption of a gas by a liquid, accompanied 
by a chemical reaction in the liquid. 

The title of Prof. Wicke’s paper 
“ Einfluss des Stofftransportes auf den Verlauf 
The discussion of 


was: 


heterogener Gasreaktionen.” 
the problem was restricted to isothermal con- 
ditions. 

First of all the concept of “Verweilzeitprozesse”’ 
They are processes in which the 
conversion, at concentration in the 
entering gas stream, is only dependent on the 
To them belong all 


was defined. 
constant 


time. 
homogeneous isothermic processes in an ideal 
ideal mixer. 
Heterogeneous processes taking place in these 


average residence 


piston-flow reactor or in an 


reactors need not be ‘ Verweilzeitprozesse,”’ 
which is the case, for instance, when a transport 
phenomenon occurs that is dependent on the 
gas velocity. 

Another definition—given by C. WaGner 
describes the ** Abklinglinge,” i.e. the length of 
a tubular reactor in which a fraction l/e of the 
reacting substance is converted. It is related to 
the concept of “height of a transfer unit”. 
are given for 
empty tubes as well as for packed beds at laminar 
and turbulent flow both when transport to the 


grains and when diffusion in the pores is deter- 


Equations for the * Abklinglinge ” 


minative, 
The velocity of a process is generally determined 
by three factors : 


(a) transport to the grains, (b) diffusion in the 
pores, (c) the reaction proper. 


As an example the oxidation of carbon with air 
in a tube with porous carbon walls was discussed 
in detail. At temperatures above 1000°C it is 
the transport in the gas to the tube wall, between 
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650 and 1000°C the diffusion in the pores and 
below 650°C the chemical reaction rate that is 
determinative. 

In the discussion Dr. Ir. A. KLINKENBERG 
objected to the definition of ‘ Verweilzeit- 
prozesse,” because actually the process is defined 
by its change. Prof. Wicker did not agree. 

Prof. Ir. H. Kramers asked why Prof. Wicker, 
in his equations, invariably uses the effective 
thickness of the boundary layer. He _ himself 
prefers transfer coefficients. Prof. WickE con- 
sidered the preference for either formulation a 
matter of taste. The use of the effective thickness 
of the boundary layer, moreover, offers certain 
advantages from a pedagogical viewpoint. 

The other questions asked during the discussion 
referred to details. 

Dr. vAN DE Vusse’s paper dealt with : 
“Engineering aspects of the oxidation of 
mercaptans in caustic solutions. Reaction kinetics 
and design of reactor.”’ In the process considered 
the liquid absorbs oxygen from the gas phase, 


after which the oxidation proceeds in the 


liquid. 


In laboratory experiments it was possible, by 
applying a high gas velocity and a high stirring 
speed, to obtain such a rapid transport of oxygen 
that the rate 
These kinetic data are used in technical reactor 


reaction proper was measured. 
design. 

First the reasons for choosing a given type of 
reactor are given. Under the prevailing conditions 
the reaction is zero-order for mercaptan and 
first-order for oxygen. It is, however, easy to 
supply air in excess. Therefore, the process need 
not be effected in counter-current and a one-stage 
stirred reactor will do. 

The optimum size of the reactor is considered 
to be the volume at which the required stirrer 
energy is minimum. OQOn_ the 
increased the reaction may proceed at a lower 
concentration of oxygen in the liquid, so that a 
lower coefficient of mass transport is required. 
At a given certain excess of oxygen, however, the 
“ superficial gas velocity” also drops when 
volume is increased, which means that the 
stirrer speed must be raised to attain the transfer 
coeflicient desired. The magnitude of the transfer 


volume being 


coefficient under given conditions of stirring is 
taken from the literature. 

The calculations lead to an optimum reactor 
volume, at which the stirrer energy is minimum. 
The concentration of dissolved oxygen is then 
about one-third of the saturation concentra- 
tion. 

In the Dr. KLINKENBERG drew 
attention to the agreement between the equations 
for absorption attended by chemical reaction 
(Harra) and the equations for pore diffusion 
(Wicke). 

Ir. Nigstnc objected to the rate of physical 


discussion 


mass transport being calculated from the measur- 
ing results of Cooper, FerNstrom and MILLER. 
These investigators measured the absorption of 
oxygen in sulphite solutions, hence, again an 
absorption with chemical reaction. In Dr. VAN DE 
Vusse’s opinion, however, in this case the physical 
rate of transport was measured all the same. 

Drs. Tu. N. ZwiererRtnG remarked that in the 
optimum situation described by Dr. van DE 
Vusse the gas velocity is very low. From an 
economical viewpoint it is probably more profit- 
able to operate at a higher gas velocity, which 
enables the stirrer energy to be reduced con- 
siderably. 

Several discussers asked for further information 
on numerical quantities used and on the deriva- 
tion of the formulae. 

The paper of Ir. Niwsinc and Prof. Kramers 
entitled : “ Absorption of CO, in carbonate — 
bicarbonate buffer solutions in a_ wetted-wall 
was presented by Ir. Nusine. 


, 


column ’ 

Experiments had been made under conditions 
at which the penetration theory is valid. This 
means that the substance to be absorbed pene- 
trates into the liquid at the same rate as with 
diffusion into a semi-infinite medium, during a 
time which is equal to the quotient of column 
height and superficial liquid velocity. As the 
concentrations of and_ bicarbonate 
ions ren.ain practically constant the concentration 
of OH ions is also constant. Thus, it is allowed 
to apply DancKWeErts’ equations to the combina- 
tion of diffusion and a first-order reaction. 

The extent to which the conditions made were 
satisfied in the experiments was ascertained in 


‘arbonate 
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detail. The only important correction taken into 
consideration was an “ end effect”. 

In applying the theoretical relation for the 
rate of absorption to the experimental data two 
unknown quantities must be taken into account, 
viz. (a) the solubility of CO, in the absorption 
liquid, which determines the CO, concentration 
at the liquid surface, and (b) the reaction rate. 
Both quantities can be deduced from the experi- 
mental data. They are both functions of the liquid 
There with 
literature data. 

In the discussion the concept of * penetration 
The 


given 


composition. is good agreement 


theory ” was brought up. definition of 
Nusinc and Kramers 


Dr. VAN DE Vusse’s opinion, however, this name 


was above. In 
is better reserved for the picture suggested by 
Danckwerts, in which turbulence “penetrates” 
as far as the phase boundary. 


Tuirp Session : “* NON-UNIFORM 


CONCENTRATION DISTRIBUTIONS ” 


Prof. Ir. H. Kramers. 


Rapporteur Tr. Ta. W. Wiuiemse. 


Chairman 


Speakers Prof. Dr. P. V. Danckwernrts (Imperial College 
of Science and Technology, London) ; Dr. K. 
Rrerema (Koninklijke /Shell 
Amsterdam); Dr. H. 
fir Chemische 
Hochschule, Darmstadt). 


Laboratorium, 
(Institut 
Technische 


HormMaNn 
Technologie, 


This session dealt with reactions taking place 
in heterogeneous or seemingly homogeneous 
systems, in which the interaction between non- 
homogeneity on the one hand and conversion and 
yield on the other play an important part. 

In Prof. Danckwerts’ paper : * The effect of 
incomplete mixing on homogeneous reactions ” 
the flow behaviour of the reactants in a continuous 
reactor For a 
first-order or pseudo-first-order reaction it is 


was discussed. homogeneous 
sufficient to know the residence-time distribution 
to enable the initial composition to be calculated, 
because the 
exclusively dependent on the duration of its 


chance of a molecule reacting is 


stay in the reactor. In the case of reactions 
with more complex kinetics the chance of a 


molecule reacting depends on the number of 


molecules it encounters during its stay, so that a 
great deal more information is required. Complete 
calculation is then only possible in a few special 
cases, e.g. the tubular reactor with longitudinal 
diffusion. 

A method of approach to the general problem 
was suggested. Two extreme cases that can 
be calculated from the residence-time distribution 


can be distinguished : 


(a) The incoming liquid is broken up into 
discrete fragments which are small com- 
pared to the reactor and in which molecules 
entering together remain together : com- 

Model : ideal tubular 

parallel, 


plete segregation. 
reactors arranged in 
case: the 
without molecular diffusion. 


The the 
molecular scale in a time which is small 


Special 


** well-stirred > tank reactor 


reactants are dispersed on 
compared to the average residence time : 
the 
mixing on the molecular scale. 


perfect mixer with homogeneous 


If the reaction mixture consists of two reactant 
streams the average reaction rate will, in general, 
depend on the ratio between these two extremes, 
‘degree of segregation”. It was 
that the 
scale between elements of fluid having different 


hence, on the 


demonstrated mixing on molecular 
ages reduces the conversion if the order of the 
reaction is higher than 1 and causes it to increase 
when this order is lower than 1. 

A quantitative definition of the degree of 
segregation was given for a two-component 
system and for the elements of different ages of 
a liquid stream passing a reactor. The requisite 
specification of the size of the areas of inhomo- 


. 


geneity was obtained by defining a “ seale of 


segregation " with the aid of a technique known 
from the statistical turbulence theory. 
The degree of segregation can be measured with 


the aid of slow reactions. In actual practice the 


study of segregation phenomena will probably 


be most useful for very fast reactions, such as 
combustion or precipitation, A few examples 
were given. 

In the discussion it was pointed out that in 
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many cases the conversion can be calculated, 
with an accuracy suflicient for practice, from the 
distribution on the 
complete segregation, viz. for slow reactions 


residence-time basis of 
of an order between } and 2 if the residence-time 
distribution more or less corresponds to that of 
a cascade consisting of more than one ideal 
mixer. 

On the other hand, with many slow reactions 
occurring in stirred vessels it will be easy to 
reach such a degree of mixing that one is sure to 
obtain homogeneity. 


A reaction is “ slow” in the sense meant here 
and a study of segregation is thus of little impor- 
the for the 


between the substance in the areas of inhomo- 


tance if relaxation time reaction 
geneity and the surrounding liquid is much 
with than the 


relaxation time for the dispersion of these areas 


longer homogeneous mixing 
by diffusion. 

Dr. Rierema entitled : 
* Heterogeneous reactions in the liquid phase. 
distribution and 


presented a paper 


Influence of residence-time 
interaction in the dispersed phase.” He discussed 
the development of heterogeneous reactions in the 
liquid phase in continuous reactors for cases in 
which the average residence times of the two 
phase do not differ much because the dispersion 
is maintained by turbulence and entrainment 
occurs. In these cases it is possible to carry 
out representative experiments in batch and, 
on the strength of their results, to predict the 
course of the reaction in the continuous reactor. 


If the reaction takes place in the continuous 


phase the reaction can be regarded as a homo- 


the 
purpose. 


treated with 
graphical for the 
This can also be done when the reaction occurs 
in the dispersed phase or in the boundary layer 
around this phase, provided the exchange in the 
disperse phase is so fast that all the particles 
of this phase invariably have the 
composition, This is possible, for instance, for 


hence be 
known 


geneous one and 


methods 


same 


reactions between two liquids, the exchange 
being effected by continuous coalescence and 
redispersion of the disperse phase. 

For the opposite phase, e.g. in a solid 
system, a new graphical method was suggested 


liquid 


and worked out for a reaction in an ideal mixer, 
in which—in DaNnckwerrts’ 
continuous phase is mixed on the molecular scale 


terminology —the 


and complete segregation takes place in the 
In the latter a distribution of 
the degree of conversion is brought about which 
distribution, 


disperse phase. 


corresponds to the residence-time 
the the 
the same throughout. 
Essentially, the method consists in the graphical 
determination 


conversion in continuous phase being 


of reaction rates, at a given 
concentration in the continuous phase and at 
various concentrations of the disperse phase, 
from the results of a series of experiments with 
different concentrations of the continuous phase. 

Complications arise with liquid —liquid reactions. 
They are caused by a difference in degree of 
dispersion between batch and continuous reac- 
tions and by an exchange to an unknown degree 
in the disperse with 
DaNcKWerts’ theory it was calculated for a 


phase. In accordance 
zero-order reaction (mass transfer fully deter- 
minative) that at 80 per cent conversion through 
complete mixing of the disperse phase the reactor 
volume can be reduced by a factor of 3. 

The influence of surface-active substances and 
contaminations is so complicated that exact 
predicting of the conversion for liquid — liquid 
reactions will probably never be possible. In the 
discussion it was mainly the behaviour of liquid 
liquid dispersions and the way in which they can 
be obtained that were brought up. 

Dr. Horman, in his paper : “* Berechnung von 
Reaktoren fiir Gemischt phase-Reaktionen”’, dealt 
with the chemical reactions which proceed to 
completion as homogeneous reactions in one or 
several phases, it being supposed that these 
phases are homogeneous or, through exchange, 
‘an be considered homogeneous. 

In many technically important reactions in the 
liquid phase, such as nitrations and sulphations, 
the reactants are intersoluble to only a limited 
degree, so that the conversion takes place in two 
(or more) phases. Consequently, the conversion 
is also dependent on the solubility and the rate of 
dissolving of the reactants in the reacting phase(s). 

An analogous case occurs when in a reaction in 
a homogeneous phase e.g. with an undesirable 
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consecutive reaction the conversion is increased 
by separating the product, for instance, by 
extraction, distillation or crystallization. 

In general, the rate of conversion will be 
determined by both the mass transfer and the 
chemical reaction. By introducing the “ mixing 
yield ” (HovGEN), i.e. the ratio of the conversion 
rate actually attained to the rate of conversion at 


ideal mixing (phase equilibrium) all cases can be 


reduced to a single one in which the chemical 
this 
value ” 


reaction is fully determinative. In case, 
made of a “distribution 


(ScHOENEMANN), Le. the reciprocal value of the 


use is 


fraction of a substance in one of the phases. 
With the aid of material balances the conversion 
calculated for a number of idealized 


apparatus, in particular also for the counter- 


was 


current case, technically interesting reactions 
serving as quantitative models. 

With the example of the production of furfural 
from an aqueous xylose solution containing 
hydrochloric acid it was demonstrated that it is 
possible to predict the optimum shape and size of 
the reactor also in complicated technical cases, 
provided the reaction rate equation is known. 

In accordance with the 


freedom : ratio between phase quantities, distri- 


many degrees of 
bution coefficient, concurrent or countercurrent 
operation, etc, there are a great many possibilities 
for effecting one and the same reaction, the 
optimum choice from which can only be made 
by checking the various possibilities through 
quantitative calculation. The non-linearity of 
the equations in case the reaction is not a first- 
order one in general calls for the use of a computer. 

In the discussion a problem with a second- 
order reaction was given which is comparatively 


easy to solve. 


* Reactor Erriciency AND 
STABILITY ” 


Fourtnu SESSION : 


Chairman Dr. Ir. A. KuUInkeENBERG. 


Rapporteur Ir. W. J. VAN ASSELT. 

Prof. Dr. K. G. Denpicn (University of 
Edinburgh and Heriot-Wall College) ; Dr. C. 
VAN Heerpen (Central Laboratorium, Staats- 
mijnen in Limburg); Dr. Ir. A. P. Ore 
(Chemiebouw, Staatsmijnen in Limburg). 


Speakers 


In this session the attention devoted 
particularly to the interaction the 
temperature and its gradient in chemical reactors 
and reactor efliciency and stability. 


was 
between 


The following cases were discussed : 
Influence of temperature on the conversion 
per unit time and per unit volume with 
simple reversible reactions. 

Influence of temperature in consecutive 
and simultaneous reactions. In this case 
a different temperature may result in a 
different vield. 
Influence of reactor 
stability. 


temperature on 


Prof. entitled : 


“ Optimum temperature sequences in reactors”, 


Denpicu, in his paper 
discussed cases (a) and (b). 

In an equilibrium reaction with a positive 
effect an temperature will 
the but shift the 
equilibrium lower conversion. A 
decrease in temperature along the length of the 


heat increase in 


increase reaction rate also 


towards a 


reactor will, in this case, cause the conversion to 
increase and it is comparatively easy to calculate 
the optimum temperature profile. The influence 
of temperature on the yield of consecutive and 
simultaneous reactions was discussed in more 
detail. 

When the activation energies of main and side 
different the 
influenced in a certain direction by changing the 
temperature. With the aid of a numerical example 
it was demonstrated that in a given case the 


reactions are reaction can be 


yield of a reaction in a two-tank cascade can be 
doubled by keeping the two reactors at different 
temperatures. 

Of course these optimum values for, the 
temperatures in the two reactors are also optimum 
values from an economical viewpoint only when 
the economy is completely determined by the 
consumption of starting material and, accordingly, 
the capital and operating costs play a minor role. 

Hence, it is possible to find an optimum 
temperature for every continuous reactor in 
which there is a given concentration of reactants 
and intermediate products. It is thus to be 
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expected that in a tubular reactor, in which the 
temperature can be controlled at any point, the 
yield can possibly be increased by continuous 
variation of the temperature over the length 
of the reactor. In the discussion Dr. F. Horn 
went further into this point. 

In the above considerations on conversion and 
yield the heat effect of the reaction was left 
out of account. In his paper, Dr. van HeEERDEN 
considered mainly those exothermic reactions in 
which the reaction heat itself is wholly or largely 
used to heat the reactants up to reaction tempera- 
ture (autothermic processes). 

In a technical continuous reactor there will be, 
at every temperature, a certain conversion and 
hence a certain quantity of heat will be liberated. 
Plotting this heat against the prevailing tempera- 
ture results in an S-curve whose first part (like 
the reaction rate) practically runs according to 
an e-power function. Afterwards, however, the 
curve is bent to a constant value of the heat 
produced, which corresponds to complete con- 
version. On the other hand the heat required 
for heating the reaction mixture to the final 
temperature, plotted as a function of this tempera- 
ture, is represented by a straight line. 

For an autothermic reaction to be realized both 
relations must be satisfied at the same time. 

In order to describe the solutions of the two 
relations Dr. vAN HEERDEN uses the quantities 
T, and AT,,,, defined as follows : 


T, is the temperature at which a conversion 
1 ¢ is obtained at isothermic piston-flow 
through a reactor with known residence 
time. (At this temperature, therefore, the 
reactor length is exactly the 
* Abklinglange ”’.) 

AT,,, is the adiabatic rise in temperature of 
the reaction mixture at complete con- 
version. 


If the reaction is carried out in an ideal mixer 

the following cases can be distinguished : 

(a) 7, (starting temperature) + AT, < T,. 
In this case the curves have only one point 
of intersection, viz. at low conversion. 
This reaction, therefore, cannot be effected 
in an ideal mixer without special measures. 


(b) T, + AT,,~T,. The curves have three 
points of intersection, two of which are 
stable. After ignition of the reaction a 
stationary state establishes itself at a 
point of high conversion and at high 
temperature. 

T, + AT, > T.. 
plete conversion at a very high tempera- 
ture. 


This results in com- 


The reaction mentioned under (a) can be 
effected in a tubular reactor with heat exchange 
between inlet and outlet. Thus, 7, is increased 
to T,' and T,’ + AT, ~ T,. 

Finally, Dr. van Heerpen considered piston 
flow with axial heat conduction, characterized 
by Péclet’s number (Pé = puc,L/A). If Pé < 1 
there is hardly any difference in temperature 
between inlet and outlet and one has an approxi- 
mation to the case of an ideal mixer with three 
stationary states (two of them stable). In the 
other extreme, Pé <1, there is only one 
stationary state. 

Whereas in the subject dealt with by Dr. van 
Heerpen and Prof. Dexsicu temperature plays 
an important role in connexion with the reaction 
velocity. Dr. Ore discussed a reaction in which 
the temperature is so high that mass transfer 
determines the overall rate. In this example 
ammonia and air (8-11 per cent NH) are passed 
over a few layers of platinum wire gauze at 
about 800°C, The heat liberated during the 
oxidation corresponds to the increase in sensible 
heat in reaction products and inert gases. 

The theoretical analogy between heat transfer 
(U) and mass transfer (x) can be written by : 


2/U 


1-03 C5* PO, 


(where C, = specific heat at constant pressure 


per kmol of gas mixture, 


P = pressure of gas mixture) 


and it is easy to deduce that in this case the 
temperature of all platinum wire 
practically equal to Ty + AT. 
When, however, the gauzes have been in use 
for some time a retardation of the reaction occurs, 


gauzes is 
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to be represented by the retardation factor #, 


which from 0-8 to 1. 
1 »1 
‘. Pr. 
The temperature of the gauzes then rises in 
the direction of flow. Besides, 
very important for burner stability. For, on the 


this factor 7 is 


gas throughput being increased and if i < 1, 
U will increase more than 2’, since x is composed 
of « and the reaction rate, of which only « 
In this the 
may be blown cold. Consequently, to eliminate 


increases. case, therefore. vauzes 


the influence of the reaction rate it is necessary, 
dependent on catalyst activity, to preheat the 
gases so as to obtain stable oxidation. 

Since ¢ results from deactivation of the catalyst, 
at 7 l preheating should be carried further. 
In addition, a temperature gradient appears in the 
gauzes, which may lead to impermissible platinum 


losses. 


Firru Session: “ Reacror DeveLorment ” 


Prof. Dr. Ir. HL. 1. WareruMan 


Chairman 
Rapporteur Dr. Tr. C. Bortmouwenr. 
Prof. Dr k. (Technische 
Hochschule, Darmstadt): Dr. F. 
(Rubrcohemie Aktiengesellschaft, Oberhausen- 
Holten), also on behalf of Dr. W. Brivz 
(Th. Goldschmidt Aktiengesellschaft, 
Mannheim-Rheinau); Dr. Ir A. Kiuinken- 
(N.N De Bataafsehe Petroleum 
Maatschappij, The Hague). 


Speakers SCHOENEMANN 


SCHNURE 


BERG 


Prof. SCHOENEMANN discussed “Die Anwendung 
der Reaktionskinetik bei der Berechnung einiger 
typischer Reaktoren”’. 

In the light of a number of examples taken 
from actual practice the importance of kinetic 
calculations for the development of chemical 
processes from the laboratory stage up to large- 
scale technical application was demonstrated. It 
the high- 


pressure hydrogenation processes for coal and 


is particularly the development of 


oil that has greatly contributed towards a clearer 
A better 


basic 


formulation of the problems involved. 


understanding was gained of simple 


concepts, such as residence time distribution in 


continuous processes and its importance for the 


computation of type, shape and size of the 
reactors, 
The 


homogeneous 


homogeneous and 


the 
divergent types of continuous reactors can be 


course of pseudo- 


reactions in most widely 
calculated if the residence time distribution and 
the 


more than one phase often the requisite experi- 


reaction rate are known: for reactions in 


mental data are lacking, particularly those 
concerning mass transfer. 

As an example of a complicated chemical 
reaction the preparation of the « xplosive hexogene 
[eyclotrimethylenetrinitramine, (CH), Ny (NOg,)s] 
by nitration of hexamethylene tetramine was 


When 
experiments a quantitative picture of the various 


discussed, through eareful laboratory 
conversions playing a part in this nitration had 
been obtained it was found possible to design at 
once, without an intermediate stage, a continuous 
technical plant which worked satisfactorily. 

In exothermic high-pressure processes often a 
circulation system is used to remove the reaction 
heat. Admittedly, this system leads to a much 
lower production per unit of volume and time, 
but it is necessary in order to be able to control 
the temperature in the reactor. This empirically 
developed method, which is in general use, can 
now be quantitatively computed in various 
cases, as was illustrated by a diagram for the 
synthesis of 1: 4-butyndiol from acetylene and 
formaldehyde with the aid of an acetylene- 
copper catalyst. 

In addition, the hydrolysis of wood according 
an example of a 


to Madison's method—as 


percolation process was discussed us We ll as 
the production of furfural from vegetable material 
containing pentosanes, where the yield is strongly 
influenced by condensation and resinification of 
the desired product under the influence of the acid 
used for the hydrolysis of the pentosanes. 

the conclusion 
the 


often 


came to 
kinetics in 


Prof. ScuoeNEMANN 
that 
calculation of 


application of reaction 


technical reactors has 


provided useful data. If these are correctly 
combined with other factors determinative of an 
exploitation that is economically justified, such 
as choice of material (also with regard to safety), 
heat economy, power consumption, ete, it is 
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possible to arrive at a more rational development 
and design of chemical works. 

In the discussion Dr. KiuINkeNBerG 
further the reaction 
hexamethylenetetramine nitration, 


went 
of the 
Prof. 
WATERMAN drew attention to recent investiga- 
Weser at Delf on the 


into mechanism 


tions carried out by 


characterization of catalysts with the aid of 


simple hyperbolic functions. 

Dr. 
Dr. Brérz, 
Berechnungen des Stoffaustausches gelegentlich 


SCHNUR presented, also on 


their paper: “* Verfahrensmassige 
der Absorption von Stickstoffoxyden bei der 
Gewinnung von Salpetersiure”. 

The rate of absorption of nitrous gases, which 
is important in nitric acid manufacture, was 
studied on the basis of measurements taken on 
a laboratory scale in a thin-film absorber, in a 
semi-technical packed column, and on a technical 
scale in a nitric acid factory. 
various factors on the rate of absorption was 
studied in order to arrive at an empirical relation 
for the The 
siderations were based on the two-tilm concept. 


mass transfer coeflicient k.. con- 

Laboratory-scale experiments showed that, in 
accordance with carlier observations of CHAMBERS 
and Surrwoop, the diffusion resistance of the 
liquid could be neglected in respect of that of the 
gas film. 

In the semi-technical absorption column filled 
with Raschig rings the mass transfer coeflicient 
was found to be dependent on the liquid through- 
put if the latter was low ; above a certain through- 
put the mass transfer coeflicient was independent 
of it if the packing is wetted completely. In 
accordance with the theory no influence of the 
diameter of the Raschig was found. The way 
of packing (rings placed or poured into the column) 
did have an influence, though. 

Although the the 
technical plant were of somewhat limited accuracy 


the 


measurements made in 


it was clear that mass transfer coetlicient 


dropped with decreasing concentration of both 


gas and acid, contrary to what was to be expected 
on the strength of the laboratory-scale experi- 
ments in the thin-film absorber. 

From the observations it might be concluded 
that the conversion proper of the nitrous gases 


behalf of 


The influence of 


with water vapour occurs in the gas phase, 
the the 
phase boundary. 

Horriuizer, NiJssine, 
Prof. vaN Kreve en took part in the discussion. 


preferably in immediate vicinity of 


Messrs. DEKKER and 


Among other things, the difficulties that may 
arise in the interpretation of the measuring results 
on account of mist formation were brought up. 

In his final remarks on the symposium Dr. 
KLINKENBERG tried to arrive at a systematic 
survey of the field of chemical reaction engineering 
in the light of the various papers presented. 

Formerly, in dealing with chemical engineering 
use was often made of a classification according to 
the so-called unit operations and unit processes. 
This symposium revealed a strong desire to go 
back to the fundamental concepts. They are 


as follows : 





Equilibria Rates 


Thermodynamics Fluid dynamics 
(including phase equilibria) Heat transport 


Mass transp< ort 


Chemical equilibria Reaction kinetics 
(simple reactions) 
Reaction kinetics 
reactions, 


(side con- 


secutive reactions, etc.) 





Chemical reaction engineering is not concerned 
the 


reactions should be 


with the chemistry of reaction, Those 


chemical classified together 
which (even though they are entirely different 
from a chemical viewpoint) are analogous as 
far as certain physical and physico-chemical 
criteria are concerned and which can be operated 
in the same type of reactor. 

The study of chemical reaction engineering is 
more diflicult than the study of processes in 
which chemical conversion does not occur (the 
unit operations). In addition to the fundamental 
concepts of flow, heat and mass transport, which 
moreover are often strongly interdependent in 
in this case, we have also to deal with chemical 
reaction. Chemical reaction engineering is more 


than a combination of unit operations and 


chemistry, as the Symposium Committee tried to 
express in the syllabus. 
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To simplify the treatment of chemical reaction 
engineering problems various methods can_ be 
applied : 


(a) Separation into equilibria and rates (or, 
if necessary, restriction to equilibria); (b) res- 
triction to flow and reaction kinetics (e.g. iso- 
thermic homogeneous reactions in stirred vessels) ; 
side and consecutive reactions may be included ; 
(c) disregard of the reaction rate in respect of the 
rate of transfer, or conversely ; (d) use of statis- 
tical considerations to simplify the treatment of 


steps that repeat themselves ; (e) use of differen- 


tial reactors. 


This was explained by referring to the various 
communications made at the Symposium. 

The theory of uniformity, too, is an important 
tool; it tries to find new solutions on the basis 
of other solutions already known. 

However, in classifying and subdividing the 
subjects belonging to chemical reaction engineer- 
ing oversimplification in formulating the problems 
should be avoided. 

Dr. Brerscunermer spoke a few words on 
behalf of the European Federation for Chemical 
Engineering to thank the Symposium Committee 


and particularly its Chairman, Dr. Ir. J. C. 
Viverer, for the considerable amount of impor- 
tant work that had been done to enable this 
twelfth meeting of the Federation. 

The closing words were spoken by Dr. D. J. 
GERRITSEN as Chairman of the Chemical Engineer- 
ing Section of the Royal Institute of Engineers. 
Dr. Gerritsen, on behalf of the two organizing 
societies, also tendered his thanks to all those 
that had taken part in the organization. He 
stressed that symposia of this nature are not only 
scientific meetings but also manifestations of a 
European collaboration which we all need so 
much nowadays. In this connexion Dr. GerrirsEN 
recalled the words spoken by Prof. Lerorr at the 
reception by the Burgomaster of Amsterdam that 
the Netherlands have always stood in the breach 
for freedom of thought and speech. 

Dr. Gerritsen then declared the Symposium 
closed. 

On the evening of the first day the participants 
were received by the Burgomaster of Amsterdam 
in the stylish surroundings of the Municipal 
Museum. On the last day the organizing societies 
invited them to a cocktail party, where the lively 
exchange of ideas was continued. 
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Notes on the diffusion-type model for the longitu- 
dinal mixing in flow (O.Lievenspiel and W.K.Smith)* 


(Received 17 September 1957) 


In their article (referred to below as “ article”) the 
authors use the diffusion type approach to characterize 
the longitudinal mixing of fluids in flow. Two methods 
are shown for evaluating the diffusion coefficient D. 
In the first method either the maximum concentration or 
the concentration at v/V = 1 (notation article and here) 
is measured. In the second method the concept of the 
variance of the concentration distribution o? is introduced. 
The article is restricted to the treatment of the infinite 
pipe: one value of D applies to the entire length and 
cross-section. 

It is the object of this letter to treat the diffusional type 
of flow in a more general manner so as to include the case 
of finite pipe length and that of varying D, applying 
suitable boundary conditions. For a discussion of the 
boundary conditions chosen we refer to an article by 
Wenner and WILHELM [1]. 

It will be shown that for the case of small D/ULT the 
influence of the special end effects can be ignored. The 
question whether the description of the mixing phenomena 
by using the mathematics of diffusion in the case of large 
D is still permissible will not be discussed. 

Futhermore it is shown that for the case where the 
average flow velocity u is known it is in some cases not 
necessary to compute the variance, as the computation 
of the mean time of passage of the tracer material leads 
to an evaluation of the diffusion coefficient D. 

In the solution use is made of the Laplace transform. 
This has the advantage of transforming the partial 
differential diffusion equation into an ordinary linear 
differential equation, the solution of which is simple. 
To derive expressions for the variance, etc. it is not 
necessary to perform the (cumbersome) back trans- 
formation. 


The differential equations 


A finite pipe of length L is bounded at both sides by a 
fore- and after-section (sections a and b resp.). As with 
Weaver and W1.HE.oM [1] these two end sections stretch 
from — oc to + . In the fore- and after-section different 
values for the diffusion coefficients compared with the 
value in the finite tube (D) occur. These are designated 
by D,, Dy resp. (Fig. 1). Injection of tracer material 
at a distance 2) downstream of the entrance of the finite 


* Article published in Chem. Engng. Sci. 1957 6 227. 


tube is supposed to take place so rapidly, without inter- 
fering with general flow, that this introduction into the 
flow can be mathematically described by the Dirac 
$-function. 

Introduction of reduced quantities and variables (see 
notation) leads to the following set of differential equations 


(la) 
=z52% (1b) 


d d 
(243-2) ane is 
38 =z) = Pep d2° 


diffusion coeff. D, D D, 
Fic. 1. 
The function R equals CV /Q. It is a kind of dimension- 


less residence time distribution. For locations downstream 
of the injection point it can be proved that (equation (10)) 


a 
[ Rae =1 22% 
0 


Boundary conditions 


In order to solve equations (la - 1c), apart from the 
obvious initial conditions <0; R, = R, = R=0, 
six boundary conditions are needed. Two conditions 
follow immediately from the consideration that at z = + « 
R, and R, must remain finite (3a, 3b). 

From the law of conservation of material at the 
boundaries there ensue two other equations (3c, 3d). 

The other two boundary equations (3e, 3f) arise from 
the intuitive argument that between the sections at the 
pipe entrance and exit the concentrations should be con- 
tinous (see also [1]). 

So the complete set of boundary equations is : 


R,(— &) = finite (3a) 


+ Here L is used as the length of the experimental section as contrasted with length of the finite tube in the rest of 


letter. 
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finite (3b) 


1 dR,(O) ; dRi(0*) 


Pe, dz Pe dz 


(3e) 


1 dR (z,>) 1 dR (z") 
° : Pe, dz 


R (2, —)— (3d) 


Pe dz 


R, (0) = R(O*) (Se) 


R (z,>) = Ry (%*) (3f) 
Properties of the Laplace transform of R 

It is convenient to introduce the Laplace transform 
of R in the form 


» 


e Pd R(d) dd (4) 


0 


R (p) 


Introduction of R stems from the consideration that 
Laplace transformation of the partial differential equations 
(la — Ie) and solving the ensuing linear ordinary differen- 
tial equations, taking into account the boundary conditions 
(3a —> 3f), leads to the solution of R as a function of z and 
p. To obtain the value for the o* it is not 
necessary to perform the cumbersome back transforma- 
tion of R(p). This can be visualized by taking the first 
two derivatives of R with respect to p and going to the 
> 0: 


variance 


limit p 
iR “ 
- | o Rad 
dp po 4 

ln | #2 Rd — b% =o? (6) 


pV ri) 


dR ) 


dp* »_0 | dp 


Equation (6), taking into account equation (5) is similar 
to equation (10) of article. 


General solution of R (p) 

Taking the Laplace transform of the partial differential 
equations (la) — (1c), solving the ordinary linear differen- 
tial equations in # (p, z), taking into account the boundary 
conditions (3a) — (3f) gives 


R,, = Aeda* 3) Pea 0 (7a) 


A nm m 
R [(q } Ga) ef 4) Pez i (q 4) Pez) : 


2q 
U 


(q = qa) € 


(= — %) [eld 4) Pe(z-z9) _. ¢-@ Pel: *0)] 


4 


R, -_ Be (db )Pe_tz =) 


where A and B are : 
A =* Pex 
4 +) areereo) 


(¢ + 4a) @ + %) A ™ — @ — ay) 


+ (q —q) e's 10) 
gPe (8a) 
~ dp) € tee 


B= R(z>2) (8b) 


The solution for z,, (= < z,, 2) is the one required : 


e (4-4) Pe 5 
R (p, Za) x 

“4 
x + Ga) F+W) + 4 


- Ga) (I — WW) 


c 2qPe (z 1 )) 4 (q Ws) € 2 Pex 


Gq) 4 ? 
+ %) — @ - Gq) (4. aqPer| 


(q + 4%q) 4 %) € 


ge 2aPe (2)-2m 





, @ — 4) 4 + YW) » 
3 
It is interesting to calculate R (p +0) as this gives the 
fraction of injected material which passed the location 2. 
From (4) and (7a —> 8b) we find (2) 
= 
R (p +0) Rd} =1 


forz = 2 


e Pe (2) for 0 


e (Pex Peg?) for = ° 


Calculation of & and o* (Table 1) 


Applications of (5) and (6) to (9) leads in a straight 
forward way to values for # and o*. The general solution 
(case I) and some special cases are listed in Table 1, 
Note that for large values of Pe the general solutions of 
% and o* tend to the following values : 


&—1 for Pe -x (lla) 


o* -»*/Pe for Pe > x (11b) 


Equation (11b) corresponds to equation (14) of article. 

The special cases were derived from the following con- 

siderations 

take values of b section equal to 
finite tube or Pe, = Pe 

>x = 1 take Pe, Pe 

© Peg = Pe = Pe, 


no diffusion in sections a and b Peg = Pe, = « 


tube froma = 0 — » 


tube from L 


tube from — ~ — 4 


The special case treated in article is case VII of Table 1. 
The variance o* in this case equals equation (12) from 
article. 


Calculations of 3 and a* from the F-curve 
Article introduced the F-curve as the output for a step 
function. It can be shown that this output is related to 
the R-curve in the following way (see also [2)}). 
9 
F(d)= | Rd db 
0 
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Letters to the Editors 


or making use of a well known relation 
transforms : 


of Laplace 


F = R/p (12a) 


Introduction of the function 


f (8) =1 — F(#) (13) 


and 


(13a) 


into equation (12a) with subsequent application of (5) 
and (6) leads to the relations 
v= }( F(d))dé 
0 
o? = 2] #(1 — F(8)) dé 


0 


Determination of the diffusion coefficient 


In those cases where D, = D, = 0 or D the equations of 
Table 1 may be used for the evaluation of the diffusion 
coeflicient D. For small values of the Peclet number Pe it 
is sufficient in many cases to have experimental information 
on & only. From the definition of # and o* there follows 
with V/o = (a, 
fluid in the experimental section, and C some concentration 


%)/u the mean residence time of the 


of the tracer, which in this case need not be dimensionless : 


x 
| aCd@ 


= 
| Cd 0 
0 

= _ Ty 
| 0? Cd 0 \@ Cd 0 \ 


0 0 \ 


o? (V/v)? —— | 8 


x ba 
| Cd @ | Cd 0 
d d 
The right hand side can be calulated from the measured 
C-curve. Furthermore $ and ¢? are known functions of the 


Peclet number (Table 1). Two cases may arise : 


Case (a): V/v is known, Now it is possible to calculate 
} immediately from (16): & =< 8/(V/v). From Table I 
8 is a known function of Pe, which value can be cal- 
culated. Calculation of o* by way of Table 1 and com- 
paring this value with s,? (v/V)* can be used as a check, 


Case (b): V/o is unknown. The calculated values 
@ and s,” are both functions of V/v and Pe. From these 
two equations V'/v and Pe can be obtained. 

To illustrate both cases, take example 2 of article : 


The value for @ and s,* are resp. : 


O = 4252 /386 sec 
8," == 48-07 sec" 


11-02 see 


Case (a): V/v = 7-7 sec. From Table 1, case VII 


b= OV v) (1 -t- - 11-02/7-7 1-452 
' Pe 


Or 1/Pe = 0-216 (article 0-217) 


Case (b) : 


- 
s,* 48-07 


(11-02)? (1 + 2/Pe)* 
1/ Pe* (2Pe + 8) 


This leads to 
0-218 (article 0-217) 


and @ | 1 1-436 makes V/o: 


» = 0/8 


== 7-67 sec (article 7-7) 


EK. Tu. van per LAAN 
N.V. de Bataafsche Petroleum Maatschappij. 
(Royal Dutch Shell Group). 


NOTATION 
Cross-sectional area of tube (L*) 
Concentration of tracer in fluid (dimensionless) 
Longitudinal diffusion coeflicient (L*/T) 
4 
t 


Rd output for a_ step function U 

§ (dimensionless) 

Length of finite tube (L) 

Parameter Laplace transformation (dimension- 

less) 

Peclet number u(xx 
= 4/(p/Pe +1 

Volume of tracer of unit concentration which 


%)  D (dimensionless) 
4) (dimensionless) 


would correspond to the actual amount of 
tracer introduced (1°) 

CV,/Q normalized concentration of tracer in 
fluid (dimensionless) 

U(%, — %)/V average flow velocity (L/T) 
Step function (dimensionless) 

Volumetric fluid flow rate (L°/T) 

A (tq 
(L*) 
Distance from entrance of finite tube (L) 
Length of experimental section of finite tube (L) 
2/(x,, — %) reduced distance from entrance of 
finite tube (dimensionless) 

D,/D ratio of diffusion coefficients in 
section and finite tube (dimensionless) 
D,/D ratio of diffusion coefficients in after 
section and finite tube (dimensionless) 

Dirac delta function (dimension @-') 

Time, measured from moment of injection (7') 


-2) volume of experimental section 


fore 
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v0/V = Ou/(x,, —% )reduced time (dimen- Subscripts 
Pertaining to fore section 


sionless) 
: Pertaining to after section 


| 3 Rd } mean of the R-curve (dimensionless) The value at point of injection 
The value at point of measurement 
The value at the end of the finite tube 


8° Rd § — 6? variance of the R curve _ bar over function designates the Laplace transform : 
o (dimensionless) = 
&(V/v) Q €PIQ(d) ds 


o* (V/v)? 
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Book Reviews 


Plant Process Dynamic Characteristics. The Pro- 
ceedings of a Conference held at Cambridge 4-6 April, 
1956 by The Society of Instrument Technology. _Butter- 
worths Scientific Publications, London, 1957. xii 
246 pp. 50s. 


Tue President of the Society, Mr. A. J. Younc, remarks 
in the foreword to this volume that, in view of the theore- 
tical and experimental work carried out in the last 10 
years, the Control Section of the Society thought it would 
be profitable to hold a conference. This book is a record of 
11 papers read with the discussions on each paper, the 
whole being welded together by a joint introduction by 
Drs. G. L. D’Omprary and G, DS, MAcLELLAN. 

In this introduction they draw attention to the advances 
made in the subject since Ivanorr’s paper to the Institute 
of Fuel in 1934 and point out that there are two classes of 
engineers interested in the subject, ic. process engineers 
who have not read or become familiar with the procedures 
common to servo-engineers and vice versa. They have 
accordingly tried to bring the two nearer together and 
perhaps it might be said that this was one of the chief aims 
of the conference. However they admit that the proceed- 
ings were weighted in the direction of chemical engineering. 

Turning now to the papers which were read, one by 
A. Horton, detailed experimental work on heat transfer in 
a jacketed pan to get information for consideration of the 
controllability of a plant in the design stage. R. L. Day 
was concerned in a similar way with the frequency response 
of a distillation column reboiler in which water was boiled 
under natural circulation in a tube 6 ft long by 1 in. in 
diameter. W. L. Wuitxinson and W. D. Armsrrone 
investigated the transient response in a 5 plate single 


bubble column separating benzene and carbon 
tetrachloride. H. Vorrrer used a laboratory Oldershaw 
sieve tray column with sections of 9, 50 and 74 trays on 
n-heptane-benzene mixtures ; he found that concentration 


formulae 


cap 


responses could be represented by simple 
containing residence times and the numbers of trays as 
the main parameters. M. W. Greer.ines investigated the 
dynamic behaviour of glass electrodes for studying pH 
in the neutralization of either 0-05 N hydrochloric or 
phosphoric acid by either 5 N or 2-5 N sodium hydroxide. 

G. H. Tuornpike, R. T. D. Ricnarnps and P. D. Jowrerr 
applied control to a woollen carding machine which reduced 
the entanglement of fibres producing a filmy web whose 
thickness was controlled by absorption of a light beam 
on a photoelectric cell. 

N. Ream, R. H. Tizanp and D. S. Townenp applied 
the N.P.L. 
cascade control loop on the heavy gas oil offtake from a 
Kent Oil Refinery crude oil still. 

J. Enxprz, J. M. VERMEULEN 
described the testing of plant units by the use of a 
pneumatic sine wave generator connected to a pneumatic 
control valve. C. A. Laws dealt with equipment for data 
reduction in which processes of automatic data logging and 


electronic analogue computer to simulate the 


L. Janssen and J. C. 


scanning were so linked that the logging equipment only 
operated when deviation was detected. Finally the volume 
includes papers which were read by J. Lewis and R. M. 
MacMinian and by J. H. Wesrcorr on mathematical 
aspects. 

Enough has no doubt been said to convince all chemical 
engineers interested in automatic control that this is a 


book which they must read. 
M. B. DonaLp 
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VOLUME 3, No. 3 SEPTEMBER, 1957 
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Tuomas J. Hanrarry and James M. ENGEN : Interaction between a turbulent air stream and a moving water surface. 
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S. K. Frreptanper : Behavior of suspended particles in a turbulent fluid. 

C.D. HoLtanp and P. G. Murpocu : Determination of the chemical properties of an oxide catalyst in a closed system. 
J. G. Stren and James H. Weser: Evaluation of the heat requirements in batch distillation operations. 

Evucene Mitier, S. P. Foster, R. W. Ross and Kurr Wont: Gas mixing in a square duct. 

N. T. Hsu and B. H. Sace: Thermal and material transfer in turbulent gas streams : local transport from spheres. 


D. L. Jounson, H. Sarro, J. D. Poteses and O. A. Hovucen: Effects of bubbling and stirring on mass transfer 
coefficients in liquids. 


J. A. Consieiio and C. M. Suiepcevicu : Effect of liquid physical properties and flow rates on the surface area of 
sprays from a pressure atomizer. 


GeorGe Tuopos: Critical constants of the aromatic hydrocarbons. 
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Belgian Technical and Chemical Engineering Journals 





J. Farcuames : Recherches sur l'enthalpie moyenne d'un écoulement gazeux. Revue Universelle des Mines, Licge 
1957 13 No. 1, 3. 

La détermination de lenthalpie moyenne d'un écoulement gazeux comporte deux mesures, celle du débit en 
poi’'s et celle de la température moyenne. Une définition rigoureuse de la température moyenne fait ressortir la 
nécessité d'une étude théorique et expérimentale des champs de vitesse et de température. 

Insuflisance d'une mesure unique au moyen des sondes classiques. Utilité d'un dispositif stabilisateur de 
lécoulement : son efficacité. Etude expérimentale du positionnement des sondes dans les écoulements symeé- 
triques et modérément dissymétriques. Confrontation des résultats avec les lois généralement admises dans 
les problémes d’écoulement et de transmission de chaleur en conduite. 


J. Vipat: L’analyse des champs thermiques par analogie électrique. Revue Universelle des Mines, Liege 1957 
13 No. 1, 32. 

Etude des problémes de conduction thermique par analogie électrique. Utilisation de la curve rhéoélectrique 
et des réseaux de résistance. Méthodes d’étude des phénoménes transitoires. Reéalisations du laboratoire de 


l'Université de Liége. 


E. Bucuer: La fluidisation. Revue Universelle des Mines, Li¢ge 1957 13 No. 3, 101. 
Etude bibliographique d’ensemble comprenant 3 parties : 

(1) rappell des données principales du calcul des réacteurs. 

(2) description des différents types d’appareils. 

(3) quelques exemples d’application industrielle. 
F. G. War.srorck : Conductibilité thermique et réactions chimiques. Bulletin de la classe des sciences d 
0’ Académie Royale de Belgique 1957 43 No. 9, 669. 
Etude théorique de l'effect de existence d'une réaction chimique dans un gaz, sur sa conductibilité thermique : 
~ influence de la pression, d'un effect catalytique de la paroi et de la forme de la cellule de mesure. 
F. Coussemant : Le Génie Chimique classique peut-il répondre aux exigences de l'industrie? Revue de la St 
Royale Belge des Ingénieurs et Industriels 1957 No. 7, 62. 


Nécessité d’étendre les méthodes du Génie Chimique aux calculs des réacteurs chimiques. 


G. de Henavu : L/industrie de raflinage du pétrole et la régulation automatique. Revue de la Sté Royale Belge 
des Ingénieurs et Industriels 1057 No. 4, 155. 


Vue d’ensemble sur les développements récents de la régulation automatique dans les raflineries de pétrole. 


J. Perers: Een vereenvoudigde metode voor de wiskundige behandeling van een geregeld industricel Process, 
aan storing onderworpen. Technisch-Wetenschappelijk Tydschrift 1957 26 15. 

Etude de la possibilité d’appliquer la méthode des diagrammes d’amplification et de phase a la régulation des 
processus industriels. 

L. Derrer, L. Liavine et P. Hesrermans : La mesure des relations P.V.T. par la méthode des isothermes. Industrie 
chimique Belge 1957 22 No. 5, 514. 

Exposé des principes et mise au point expérimentale d'un dispositif d’étude de la compressibilité des gaz industriels 
entre 0 et 200°C et jusque 3000 atm d’aprés la méthode d’Amagat. 

J. Caruara : L’ingénieur du Génie Chimique et les besoins de l Industrie. Industrie chimique Belge 1957 22 No. 6, 
665. 

Conférence inaugurale dun cycle de cours post Universitaires. Définition exacte du Génie Chimique — Son 
intérét pratique. Neécessité de lorganisation de cours post universitaires. 

M. L. Hetumeckx: Echange de matiére entre phases. Industrie chimique Belge 22 No. 8, 907 1957. 


Principes des méthodes actuelles d’étude des processus d’échange de matiére entre phases. 


A. Peyravy et P. Focu : Etude de la semi carbonisation par fluidisation. Résultats obtenus dans une installation 
pilote. Annales des mines de Belgique 1957 No. 1, 32. 
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31 1999. 
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ERRATA 


Hoetscner H. E., Temperature stability of fixed-bed catalytic converters. Chem. Engng. Sci. 1957 6@ 183. 
Equation (22) in top right-hand column of p. 186 should read : 


10 
T= (mo + - 


Equation (33) on p. 188 should read as follows : 


( 


As indicated in equation (5) : 


For case ITI 


where f{(p) = 


Thus y 


Hamer, Eom and Lypersen A. L., The vapour pressure of di-n—butylphthalate, di-n—butylsebacate, laurie acid 
and myristic acid. Chem. Engng. Sci. 1957 7 66. 

On p. 69 Figs. 1 and 2 should be transposed. 

Barren M. H. I., McGrar W., Rumrorp F. and Siesser C. G. M., Some considerations on heat transfer in spiral plate- 


heat exchangers. Chem. Engng. Sci. 1957 7 112. 


Fig. 1, p. 112 should be marked A, B, C and D, as in the revised figure below. 
12° -—- 
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